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1.0  Background 


The  key  objective  of  the  Joint  Advanced  Distributed  Simulation  Joint  Test  Force  (JADS  JTF)  is 
to  provide  the  Test  and  Evaluation  (T&E)  community  with  an  evaluation  of  the  utility  of 
Advanced  Distributed  Simulation  (ADS)  as  a  methodology.  The  End  To  End  (ETE)  test 
evaluates  the  utility  of  ADS  to  complement  the  Developmental  Test  and  Evaluation  (DT&E)  and 
Operational  Test  and  Evaluation  (OT&E)  of  a  Command,  Control,  Communications,  Computers, 
Intelligence,  Surveillance  &  Reconnaissance  (C'lSR)  system.  The  test  uses  the  critical 
operational  issues  for  the  Joint  Surveillance  Target  Attack  Radar  System  (Joint  STARS)  to 
conduct  its  T&E  utility  evaluation  in  an  ADS -enhanced  test  environment. 

The  ETE  test  consists  of  four  phases.  The  first  two  phases  occur  in  a  laboratory  environment, 
suited  for  exploring  DT&E  and  early  OT&E  applications.  Phase  3  checks  compatibility  of  the 
ADS  environment  with  the  actual  Joint  STARS  equipment,  and  Phase  4  is  a  live  open-air  test 
designed  to  mix  live  and  virtual  targets  and  provide  an  end-to-end  environment  for  testing  Joint 
STARS  in  its  operational  environment.  The  intent  is  to  provide  a  set  of  interfaces  from  sensor  to 
weapon  system  including  some  of  the  intermediate  nodes  that  would  be  found  in  a  tactical 
engagement.  The  test  traces  a  thread  of  the  battlefield  process,  from  target  detection  to  target 
assignment,  target  engagement,  and  battle  damage  assessment  at  corps  level,  using  ADS.  It 
allows  the  tester  to  evaluate  the  thread  as  a  whole  and  to  evaluate  what  effects  an  operationally 
realistic  environment  has  on  the  system  under  test.  The  ETE  test  is  designed  to  add  additional 
entities  in  a  seamless  manner  to  the  battlefield  seen  by  Joint  STARS.  In  addition,  the  ETE  test 
adds,  via  ADS,  some  of  the  complimentary  suite  of  the  C''lSR  systems  and  weapons  systems  with 
which  Joint  STARS  interacts.  This  enables  the  test  team  to  evaluate  the  utility  of  an  ADS- 
enhanced  test  environment. 

The  ETE  test  uses  ADS,  as  defined  by  IEEE  Standard  1278  for  Distributed  Interactive  Simulation 
(DIS),  to  supplement  the  operational  environment  that  E-8C  and  Light  Ground  Station  Module 
(LGSM)  operators  would  experience.  By  mixing  any  available  live  targets  with  targets  generated 
by  a  simulation,  the  ETE  synthetic  environment  presents  a  battle  array  that  represents  many  of  the 
major  ground  systems  found  in  a  corps  area  of  interest.  Additionally,  by  constructing  a  network 
with  nodes  representing  appropriate  ClSR  systems  and  weapon  systems,  a  more  robust  cross 
section  of  players  is  available  with  which  the  E-8C  and  LGSM  operators  can  interact. 

Several  components  are  required  to  create  the  ADS-enhanced  operational  environment  (ETE 
synthetic  environment)  that  is  used  in  the  ETE  test.  In  addition  to  Joint  STARS,  the  ETE  test 
requires  Janus,  a  validated  simulation  capable  of  generating  thousands  of  entities  that  represent 
some  of  the  elements  in  a  threat  rear  area  of  operation.  Also,  simulations  of  the  Joint  STARS 
moving  target  indicator  (MTI)  radar  and  synthetic  aperture  radar  (SAR),  collectively  called  the 
Virtual  Surveillance  Target  Attack  Radar  System  (VSTARS),  are  used  to  insert  the  simulated 
entities  into  the  radar  stream  aboard  the  E-8C.  Other  components  used  to  support  the  test  include 
live  elements  of  the  Army’s  artillery  command  and  control  process  and  the  Tactical  Army  Fire 
Support  Model  (TAFSM),  which  simulates  a  Battalion  of  the  Army’s  Advanced  Tactical  Missile 
System.  These  simulations  begin  to  interact  after  an  operator  starts  a  scenario  in  the  DIS  version 
of  Janus.  As  VSTARS  processes  the  simulated  entities,  the  LGSM  receives  MTI  and  can  request 
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SAR  images.  Using  doctrinally  correct  means,  a  soldier  sends  free  text  messages  from  the 
Compartmented  All-Source  Analysis  System  Message  Processing  System  to  two  remote 
workstations  (RWSs).  In  turn,  a  soldier  at  an  Advanced  Field  Artillery  Tactical  Data  System 
(AFATDS)  receives  Target  Intelligence  Data  messages  from  the  RWSs.  The  AFATDS  operator 
sends  a  fire  mission  to  another  AFATDS  operating  as  a  Battalion  Fire  Direction  Center.  Here, 
TAFSM  encapsulates  the  fire  and  detonation  traffic  within  DIS  protocol  data  units  (PDUs)  and 
broadcasts  the  PDUs  across  the  ETE  synthetic  environment.  Finally,  Janus  receives  the  PDUs, 
assesses  damage,  and  continues  the  end  to  end  loop. 

2.0  Introduction 

With  the  rapidly  changing  environment  of  distributed  testing,  one  thing  must  stay  constant — 
digital  terrain.  In  September  1996,  the  ETE  test  team  realized  that  the  simulations  being 
developed  to  support  ADS  needed  to  work  from  the  same  terrain  database.  In  order  to  reduce 
the  level  of  effort  required  to  develop  digital  terrain  for  the  ETE  test,  the  team  used  the  relatively 
featureless  terrain  of  Southwest  Asia  (SWA). 

Three  of  the  four  simulations  (Janus,  TAFSM,  and  the  VSTARS  MTI)  in  the  ETE  synthetic 
environment  already  used  National  Imagery  and  Mapping  Agency  (NIMA)  digital  data. 
Therefore,  the  ETE  team  decided  to  use  this  terrain  data  as  a  basis.  The  team  contracted 
Lockheed  Martin  Tactical  Defense  Systems  (LMTDS)  of  Litchfield  Park,  Arizona  to  design  and 
build  the  fourth  simulation,  a  Joint  STARS  SAR  simulation  called  the  Advanced  Radar  Imaging 
Emulation  System  (ARIES). 

From  a  distributed  testing  standpoint,  it  was  essential  that  each  simulation  represent  terrain  and 
features  accurately.  Therefore,  all  four  simulations  used  terrain  data  derived  from  Level  1  Digital 
Terrain  Elevation  Data  (DTED)  and  Digital  Feature  Analysis  Data  (DFAD).  Despite  using  a 
common  terrain  basis,  approximately  90%  of  the  development  effort  was  correcting  the  data, 
adding  detail,  and  putting  it  into  a  standard  format.  This  was  accomplished  by  combining  the 
strengths  of  Geographic  Information  System  (GIS)  software  tools  such  as  ArcView®  and 
ARC/INFO®  with  the  flexibility  of  PV-WAVE®,  and  C  and  FORTRAN  code. 

The  process  and  standards  involved  in  producing  terrain  databases  needs  to  change  to  make  it  less 
time  and  manpower  intensive.  This  report  focuses  on  the  steps  used  to  develop  the  ETE  terrain 
database,  lessons  learned,  and  what  can  be  done  to  improve  the  whole  process  in  the  future. 

3.0  Methodology 

The  ETE  test  team  based  its  testing  on  a  54-hour  Corps  Battle  Simulation  scenario  in  SWA,  used 
by  TEXCOM  Lab  for  testing  C"*!  systems.  This  scenario  was  adapted  from  the  US  Army 
Command  and  General  Staff  College  (CGSC)  Common  Teaching  Scenario  -  Southwest  Asia, 
dated  April,  1992,  modified  by  Headquarters,  TRADOC.  Not  only  did  the  scenario  dictate  the 
entities  present  in  Janus,  but  it  also  defined  SWA  as  the  area  of  interest  for  the  ETE  terrain 
database.  Therefore,  in  order  to  interact  in  the  ETE  synthetic  environment,  Janus,  TAFSM,  and 
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VSTARS  required  a  terrain  database  of  SWA.  The  following  steps  highlight  the  major  concepts 
used  in  developing  the  ETE  terrain  database. 

3.0.1  Designing  the  Advanced  Radar  Imaging  and  Emulation  System 

ARIES,  a  component  of  VSTARS,  ultimately  defined  how  the  ETE  terrain  database  would  be 
built.  During  the  ARIES  design  process,  two  implementations  were  considered — a  point  map  and 
raster  ground  truth  generation.  Figure  1  represents  an  overview  of  the  final  A  RTFS  system 
design.  As  shown  below,  the  point  map  design  was  selected.  At  the  time,  this  design  offered 
distinct  advantages.  First,  DTED  and  DFAD  data  was  readily  available.  Second,  the  tools  to 
modify  and  create  new  terrain  databases  were  available.  Finally,  all  of  the  simulations  in  the  ETE 
synthetic  environment  could  derive  their  own  proprietary  formats  from  the  single  terrain  database. 
Refer  to  Appendix  C  for  the  interface  control  document  (ICD),  which  defines  the  characteristics 
and  structure  of  the  digital  terrain  elevation  and  feature  database  for  the  ARIES  simulation. 
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3.0.2  Getting  the  Data 


The  EXE  team  obtained  terrain  data  as  paper  maps  and  in  digital  format  from  NIMA.  The 
Catalog  of  Maps,  Charts,  and  Related  Products  and  the  Semiannual  Bulletin  Digest  proved  to  be 
invaluable  for  determining  availability  of  terrain  data  for  SWA.  Figure  2,  from  TEXCOM’s  The 
Road  to  War,  shows  the  area  of  interest  for  ETE  terrain  database  development. 
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3.0.2. 1  National  Imagery  and  Mapping  Agency  Digital  Data 


Level  1  DTED  and  DEAD  can  be  obtained  on  CD-ROM  from  the  National  Imagery  and  Mapping 
Agency,  3200  S.  Second  St.,  St.  Louis,  MO  63118.  The  request  should  indicate  full  one-degree 
geographic  cells  of  data  by  either  southwest  comer  coordinates,  or  by  delineating  the  required 
area. 

3.0.2.2  Digitizing 

Since  the  Level  1  DEAD  did  not  provide  sufficient  detail  for  ARIES,  TEXCOM  labs  obtained 
medium  scale  paper  maps  (e.g.,  1:50,000  and  1:100,000)  from  the  NIMA.  The  ETE  test  team 
digitized  features  from  these  maps  to  augment  Level  1  DEAD. 

3.0.3  Putting  the  Data  into  a  Standard  Form 

When  pulling  the  data  together,  it  was  unlikely  that  it  would  all  be  in  a  consistent  format. 
Therefore,  three  important  references  were  considered  during  ETE  terrain  database 
development — datums,  map  projections,  and  map  scales. 

3.0.3. 1  Datums 

In  general  terms,  the  size  and  shape  of  the  Earth  is  modeled  as  a  spheriod.  A  geodetic  datum  uses 
this  approximation  to  define  the  mathematical  relationship  of  the  size  and  shape  of  the  Earth  to  a 
coordinate  system.  Since  there  are  many  methods  to  describe  the  Earth,  different  countries  and 
agencies  use  different  datums  to  identify  coordinates  in  GIS  software.  Referencing  coordinates  to 
the  wrong  datum  can  result  in  position  errors  of  several  kilometers,  so  it  is  important  to  be  aware 
of  the  variety  of  datums.  Coordinates  in  the  ETE  terrain  database  are  in  World  Geodetic  System 
1984  (WGS84). 

3.0.3.2  Map  Projections 

A  map  projection  portrays  the  surface  of  the  Earth  on  a  2-D  plane.  Unfortunately,  projection 
always  creates  distortions  of  the  Earth’s  surface  in  shape,  scale,  and  area.  The  key  to  selecting 
the  best  projection  is  determining  which  projection  minimizes  those  distortions  most  important  to 
the  cartographer.  DEAD  is  distributed  in  geographic  coordinates  (latitude  and  longitude).  All  of 
the  paper  maps  used  in  developing  the  ETE  terrain  database  were  in  the  Universal  Transverse 
Mercator  (UTM)  projection.  In  order  to  merge  the  two  in  GIS  software,  UTM  coordinates  were 
transformed  into  geographical  coordinates.  Later,  the  coordinates  were  projected  into  the 
Topocentric  Coordinate  System  (TCS)  for  use  in  VSTARS. 

3.0.3.3  Map  Scales 

The  map  scale  (e.g.,  1:1,000,000  or  1:50,000)  is  the  relationship  between  the  distance  on  a  paper 
map  to  the  same  distance  on  the  Earth’s  surface.  Mathematically,  this  relationship  is: 
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_ 1 _  =  rMap  Distance")  x  (Units  Conversion) 

(Scale  Denominator)  (Earth  Distance) 

The  map  scale  also  determines  how  features  are  depicted  on  a  paper  map.  Feature  representation 
changes  with  map  scale.  A  map  scale  with  a  small,  scale  denominator  has  the  greatest  feature 
detail  and  is  ironically  considered  to  be  a  large-scale  map.  With  a  small-scale  map,  the  location 
and  size  and  shape  of  features  become  distorted.  Some  features  are  even  omitted.  Therefore,  it  is 
best  to  work  with  the  largest  map  scale  that  is  available. 

3.0.4  Manipulating  the  Data 

The  first  question  in  developing  a  digital  terrain  database  is  “How  large  of  an  area  is  the  database 
going  to  cover?”  This  is  an  important  question  in  that  the  answer  determines  how  the  terrain 
database  is  built.  A  second,  and  equally  important  question,  is  “What  is  the  location  of  the 
database?”  An  area  with  a  number  of  features  also  affects  how  the  terrain  database  is  built.  In 
either  case,  the  idea  is  to  convert  individual  DFAD  cells  into  the  largest  area  possible. 

A  DFAD  file  consists  of  a  set  of  manuscripts  that  contain  point,  line,  and  area  features  over  a  1° 
by  1°  geographic  cell.  Each  manuscript  is  a  database  of  geographic  coordinates  and  attributes  that 
identify  natural  and  man-made  features  to  a  specific  level  of  accuracy.  See  MIL-PRF-89005  for 
the  DFAD  performance  specification.  When  ARC/INFO®  converts  a  DFAD  file,  it  creates  a 
workspace  with  two  coverages  per  manuscript  (DS01P...DSOnP  and  DS01L...DS0nL),  where  n 
corresponds  to  the  number  of  manuscripts  in  the  DFAD  file.  See  Figure  3  for  the  graphical 
relationship  between  an  ARC/INFO®  workspace  and  DFAD  file. 


ARC/INFO  Workspace 


coverage  DS01L 
(lines  &  points) 

coverage  DS01P 
(areas) 


coverage  DS02L 
(lines  &  points) 

coverage  DS02P 
(areas) 


ho 


DFAD  File 


Manuscript  #1 


Manuscript  #2 


A - DFADARC 

ARCD  FAD - ► 


Figure  3  Converting  between  ARC/INFO®  and  DFAD 
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Each  DSO«L  coverage  contains  both  line  and  point  features.  Each  dsOmP  coverage  contains  area 
features.  Furthermore,  ARC/INFO®  coordinate  and  attribute  files  are  related  by  an  item  called 
FACODE.  The  FACODE  in  the  DSO«L.ACODE,  DSO«L.XCODE,  and  DSO«P.PCODE  files  are  related 
to  the  DSOnL.AAT,  DSO/jL.PAT,  and  DSOnP.AAT  cover_IDs  respectively.  Figure  4  illustrates  how 
point,  line,  and  area  features  and  attributes  are  stored  in  ARC/INFO®. 


DFADFile  ARC/INFO  Workspace 


Figure  4  Relationship  between  DEAD  and  ARC/TNFO®  Attributes 

The  following  paragraphs  give  a  technical  “how-to”  create  digital  terrain  databases.  Figure  5 
illustrates  the  ETE  terrain  database  development  process.  Of  particular  note,  the  “edgematch  & 
attribute  check”  and  “digitize”  steps  require  the  greatest  resources.  Although  this  is  the  method 
used  to  create  the  ETE  terrain  database,  there  may  be  a  better  process.  The  recommendations 
section  in  this  report  explores  other  techniques  that  may  be  used  in  the  future.  Since  the  terrain 
database  was  created  using  ARC/INFO®,  software  specific  command  line  entries  are  bolded  with 
square  brackets.  These  steps  also  assume  correct  references  to  datum,  projection,  and  scale. 
Finally,  the  source  code  for  any  references  to  ARC/INFO®  macro  language  (AML)  scripts,  PV- 
WAVE®  procedures,  and  C  and  FORTRAN  code  can  be  found  in  the  appendices  of  this  report. 

RunDFADARC. 


[DFADARC  <dfad  file>  <workspace>]  This  command  converts  the  DFAD  cells 
that  will  be  used  in  the  ETE  terrain  database  into  ARC/INFO®  coverages  and 
workspaces. 


Join  attributes. 


[run  join_codes.amI]  This  AML  joins  ARC/TNFO®  attribute  files  with  the 
corresponding  ARC/INFO®  feature  file.  Refer  to  Figure  4  to  see  how  files  are 
related. 

Create  large  coverage  (s). 

This  step  creates  a  large  coverage  fi-om  several  ARC/INFO®  coverages.  A  single 
coverage  is  created  for  the  area  features,  and  a  single  coverage  is  created  for  the  line  and 
point  features. 

Area 

[put  <coverage>]  This  command  is  performed  for  each  coverage  that  will 
comprise  the  large  coverage.  Select  the  boundary,  put  it  to  a  new  coverage  called 
MASTERBNDS,  and  then  delete  the  boundary. 

[append  <new  coverage>  line  features]  This  command  joins  the  area  coverages 
into  a  large  coverage. 

[BUILD  <new  coverage>  line]  This  command  builds  area  topology  for  the  new 
coverage. 

[CLIP  <  new  coverage  >  <clip  box>  <clip  coverage>  line  0.00001]  If  a  subset  of 
the  large  coverage  is  desired,  this  command  will  “cookie  cut”  the  large  coverage 
based  on  the  shape  of  the  clip  coverage.  After  using  this  command,  manually 
“close”  any  polygons  that  have  been  clipped. 

[BUILD  <clip  coverage>  line]  This  command  builds  topology  for  the  new  area 
coverage. 

Lines  &  Points 

[APPEND  <new  coverage>  link  features]  This  command  joins  the  line  and  point 
coverages  into  a  large  coverage 

[BUILD  <new  coverage>  line]  This  command  builds  line  topology  for  the  new 
coverage. 

[build  <clip  coverage>  point]  This  command  builds  point  topology  for  the  new 
coverage. 

[CLIP  <new  coverage>  <clip_box>  <clip  coverage>  link  0.00001]  If  a  subset  of 
the  large  coverage  is  desired,  this  command  will  “cookie  cut”  the  large  coverage 
based  on  the  shape  of  the  clip  coverage. 

[build  <clip  coverage>  line]  This  command  builds  line  topology  for  the  new 
coverage. 

[BUILD  <clip  coverage>  point]  This  command  builds  point  topology  for  the  new 
coverage. 
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Edgematch  &  attribute  check. 

First,  manually  “edgematch”  each  new  coverage  with  a  back  coverage  of 
MASTERBNDS.  MASTERBNDS  provides  a  template  for  all  of  the  original  cell 
boundaries.  It  is  useful  because  it  identifies  the  areas  that  need  to  be  edgematched. 

Next,  correct  or  create  attributes  for  each  feature.  See  appendix  A  for  a  list  of  supported 
features  in  ARIES.  Appendix  B  shows  the  values  used  for  attributes  of  features  in  the 
ETE  terrain  database. 

Each  feature  needs  a  distinct  ID.  The  following  commands  create  separate  IDs  for  all 
features  in  the  ETE  terrain  database. 

Use  these  commands  for  DSOnL.AAT. 

[CALC  dsOll-id  =  dsOll-id  +  10000] 

[CALC  facode  =  dsO  11-id] 

[IDEDIT  dsOll  line] 

Use  these  commands  for  DSOwL.PAT. 

[CALC  dsOll-id  =  dsOll-id  +  20000] 

[CALC  facode  =  dsOll-id] 

[IDEDIT  dsOll  point] 

Extract  the  attributes  from  each  coverage  as  temporary  .acode,  .pcode,  and  .xcode  files. 

[INFODBASE  <info  file>  <dbase  file>]  This  command  converts  ARC/INFO®  files 
into  dbase  format. 

Use  Microsoft®  Excel  to  read  the  dbase  files.  Eliminate  the  duplicate  feature  entries. 
NOTE:  The  boundary  feature  should  always  have  an  ID  =  1. 

[DBASEINFO  <dbase  file>  <info  file>  define]  This  command  converts  dbase  files 
back  to  ARC/INFO®  format. 

.acode  Define  the  ARC/INFO®  line  files  as  shown  below. 

facode  facode  4  5  B 

height  height  4  5  B 

ficode  ficode  4  5  B 

smccode  smccode  4  5  B 

direct  direct  111 

width  width  2  3  B 

.xcode  Define  the  ARC/INFO®  point  files  as  shown  below. 

facode  facode  4  5  B 

height  height  4  5  B 

ficode  ficode  4  5  B 

smccode  smccode  4  5  B 

orientatio  orientation  2  3  B 

length  length  2  3  B 
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width  width  2  3  B 

.pcode  Define  the  ARC/INFO®  area  files  as  shown  below. 

facode  facode  4  5  B 
height  height  4  5  B 
ficode  ficode  4  5  B 
smccode  smccode  4  5  B 
nstruct_ps  nstruct_psk  2  5  B 
pct_tree_c  pct_tree_cov  2  5  B 
pct_roof_c  pct_roof_cov  2  5  B 

[RUN  modify_codes.aml]  This  AML  organizes  .pcode,  .acode,  and  .xcode  files  so 
that  features  are  incremented  by  1  in  ascending  numerical  order. 


Figure  5  Terrain  Database  Development  Process 


Creating  ARC/INFO®  coverages  from  ArcView®  projects. 

This  step  gets  digitized  features  from  ArcView®.  After  this  conversion,  the  previous  step, 
edgematch  &  attribute  check,  must  be  redone. 

[SHAPEARC  <shape  file>  <coverage>]  This  command  converts  an  ArcView® 
shape  file  into  an  ARC/INFO®  coverage. 

Run  ARCDFAD. 

[arcdfad  <workspace>  <dfad  file>]  This  command  converts  ARC/TNFO® 
coverages  into  a  DFAD  file. 

Converting  from  DFAD  to  ARIES  format. 
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dfadlbits.f  This  FORTRAN  program  converts  32  bit  DFAD  to  36  bit  words  and 
writes  an  ASCII  file  of  36  bit  words. 

bits2aries.f  This  FORTRAN  program  reads  the  36  bit  ASCII  file  and  writes 
DFAD  features  to  ARIES  format  but  with  coordinates  in  latitude-longitude. 

aries2tcs.pro  This  PV-WAVE®  procedure  converts  latitude-longitude  coordinates 
into  TCS  coordinates.  The  ARIES  file  format  is  TCS  (X,  Y)  only. 

tcs2bin.pro  This  PV-WAVE®  procedure  eonverts  the  ARIES  file  to  binary 
format. 


Figure  6  Converting  to  the  ARIES  format 


Creating  Janus  terrain  from  ARC/INFO®  coverages. 

[UNGENERATE  <line  |  point>  <coverage>  <text  file>]  This  command  converts 
ARC/INFO®  coverages  into  (x,  y)  point  files. 

janus2arc.c  This  C  program  converts  (x,  y)  point  files  into  a  format  that  can  be 
read  by  Janus  terrain  input  tools. 

3.0.5  Limitations 

The  greatest  limitation  of  creating  the  ETE  terrain  database  is  the  quality  of  the  data.  Although 
most  of  the  digital  data  is  obtained  from  NIMA,  there  is  no  consistency  between  geographic  cells 
of  data.  Often,  features  shown  on  a  comparably  scaled  paper  map  do  not  match,  or  worse,  are 
missing.  This  brings  up  an  interesting  dilemma.  “Do  you  use  DFAD  Level  1  in  its  original 
format?”  In  the  case  of  the  ETE  terrain  database,  more  detailed  digital  data  was  not  available. 
“Or,  do  you  spend  time  and  resources  to  digitize  the  detail,  correct  features  across  geographic  cell 
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boundaries,  and  modify  feature  attributes  to  make  them  more  realistic?”  That  is  what  happened 
with  the  ETE  terrain  database. 

Only  16,383  features  are  allowed  in  each  DEAD  manuscript.  Since  the  ETE  terrain  database 
ultimately  contained  over  45,000  features,  a  single  manuscript  was  not  possible.  Since  each 
manuscript  contains  lines,  points,  and  area  features  (refer  to  figure  4),  it  made  sense  to  create 
several  smaller  manuscripts  to  define  the  area.  The  basic  premise  was  to  divide  manuscripts  based 
on  existing  DEAD  geographic  cell  boundaries.  This  organized  the  area  and  made  it  easier  to 
edgematch.  Dividing  cells  independently  of  existing  boundaries,  while  possible,  would  have 
added  unnecessary  complexity  and  manhours  to  the  project.  The  benefit  of  this  approach  was  that 
it  was  possible  to  set  up  separate  manuscripts  to  distinguish  between  original  DEAD  and  the 
digitized  detail. 

The  computer  hardware  (Hewlett-Packard  735  workstation  with  144-megabyte  memory  and  8- 
gigabyte  hard  disk  storage)  used  to  run  ARC/INEO®  was  adequate  for  the  project;  however,  using 
a  comparatively  inferior  workstation  could  jeopardize  a  user’s  ability  to  quickly  navigate  (zoom 
in,  out,  or  across)  coverage  views.  Computer  disk  space  was  never  a  concern.  The  entire  ETE 
terrain  database  could  have  been  completed  using  a  4-gigabyte  hard  disk  drive.  A  Pentium®  class 
personal  computer  and  digitizer  is  required  for  digitizing  into  ArcView®  and  performing  other 
data  manipulations,  finally,  the  computers  need  to  be  networked  to  facilitate  the  exchange  of 
data  between  ArcView®,  ARC/INEO®,  and  Microsoft®  Excel. 

finally,  during  the  ETE  terrain  database  development,  there  was  no  ability  to  view  features  as 
they  would  appear  in  ARIES,  feature  attributes  (length,  width,  height),  as  defined  by  the  DEAD 
specification,  were  coded  in  the  database  as  half  their  actual  value  and  then  only  as  a  whole 
number.  Also,  the  feature  orientation  was  coded  in  the  database  as  a  whole  number  from  zero  to 
eight,  with  each  number  representing  the  orientation  of  the  feature’s  length  as  a  multiple  of  11.25 
degrees  from  true  north.  Refer  to  appendices  A  and  B  for  the  following  example.  A  single  family 
dwelling  (coded  as  length  =  6,  width  =  4,  height  =  2,  and  orientation  =  4)  would  be  a  12-meter  by 
8-meter  by  4-meter  feature  oriented  at  45  degrees  in  ARIES.  Imagine  thousands  of  features  in 
any  given  area.  It  would  be  nice  to  know  how  these  attributes  affect  one  another  with  a  3-D 
visual  representation. 

4.0  Lessons  Learned 

In  order  to  accomplish  the  project,  it  is  essential  to  not  only  have  the  software,  but  have  personnel 
trained  to  use  it.  It  took  approximately  18  months  to  develop  the  ETE  terrain  database,  four  of 
those  months  were  spent  learning  the  GIS  software  and  developing  the  processes  that  would  be 
used  to  create  the  database.  An  individual  with  GIS  software  experience  could  easily  shave  25% 
off  this  timeline,  if  not  more.  Adding  additional  people  would  also  decrease  development  time. 
In  this  case,  resources  would  be  best  allocated  during  the  laborious  edgematching  and  digitizing 
steps  shown  in  figure  5. 
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The  level  of  effort  required  to  develop  a  terrain  database  is  time  and  manpower  intensive,  but 
careful  planning  can  shorten  the  development  schedule.  Time  can  be  saved  by  not  having  to 
correct  the  digital  data. 

•  Get  the  most  recent  and  accurate  data  available  from  NIMA. 

•  Use  a  consistent  datum,  for  example  WGS84,  when  digitizing  from  paper  maps. 

•  Determine  what  fidelity  is  required  for  the  terrain  database. 

The  ETE  terrain  database  was  developed  for  the  relatively  featureless  desert  of  SWA.  The 
additional  terrain  and  feature  complexity  of  Bosnia  or  Korea  would  require  a  GIS  team  and 
systematic  plan  to  divide  and  conquer  the  terrain  database. 

Map  error  is  cumulative  and  comes  from  various  sources.  Using  map  projections  introduces 
error.  Carelessly  mixing  datums  introduces  error.  Digitizing  introduces  error.  Converting 
between  file  formats  of  different  terrain  databases  introduces  error.  The  only  way  to  minimize 
map  error  in  distributed  simulations  is  to  ensure  that  all  simulations  use  a  common  terrain 
database. 

5.0  Recommendations 

There  are  several  emerging  technologies  that  could  improve  ETE  terrain  database  development  in 
the  future.  First,  an  easy  method  to  transfer  features  from  intelligence  sources  (such  as  satellite 
photos)  into  a  digital  format  would  provide  the  capability  to  quickly  add  exceptional  detail  to 
terrain  databases.  R2V"'  from  Able  Software  Company,  ERDAS®  IMAGINE  Advantage™  with 
IMAGINE  vector  module,  and  AUTOGRAPHICS®  from  LMTDS,  Akron,  OH  are  three  software 
packages  that  would  facilitate  raster  to  vector  conversion  during  terrain  database  development. 
Second,  eliminating  conversions  between  file  formats  of  different  terrain  databases  would 
minimize  map  errors  among  distributed  simulations.  Vector  Product  Format  (VPF)  is  meta-data 
and  the  Synthetic  Environment  Data  Representation  &  Interchange  Specification  (SEDRIS)  is  a 
meta-model  that  promises  to  be  compatible  with  a  wide  variety  of  applications. 

R2V™  is  a  simple,  intuitive  software  package  that  automatically  vectorizes  raster  images.  R2V™ 
has  standard  vector  editing  tools,  but  its  sub-par  image  processing  tools  would  limit  projects  with 
complex  images.  R2V™  can  label,  georeference,  and  export  vector  data  to  other  major  GIS 
formats.  Another  great  feature  is  the  ability  to  merge  multiple  vector  files.  For  more  information, 
visit  Able  Software  Company  on  the  web  at  http://www.ablesw.com. 

The  ERDAS  IMAGINE®  software  package  is  more  robust  than  R2V™.  ERDAS  IMAGINE® 
works  with  a  variety  of  raster  and  vector  formats.  ERDAS®  IMAGINE  Advantage™  allows  the 
user  to  directly  access  image  data  in  native  format,  and  then  display  and  link  multiple  data  files. 
When  combined  with  the  vector  module,  this  important  feature  gives  the  user  the  ability  to  create 
and  edit  vector  data  from  images.  Orthocorrection,  advanced  image  processing,  and  spatial 
analysis  make  it  easier  to  develop  accurate  databases.  For  more  information,  visit  ERDAS®  on 
the  web  at  http://www.erdas.com. 
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The  AUTOGRAPHICS®  software  package  allows  the  user  to  “train”  the  software  to  extract 
features  fi’om  a  raster  image  of  a  paper  map.  First,  the  user  selects  an  example  of  each  feature 
from  the  raster  image  with  point-and-click  actions.  After  the  user  identifies  examples,  the 
software  automatically  classifies  the  remaining  features  in  the  image.  For  more  information, 
contact  LMTDS  Business  Development  in  Akron,  Ohio  at  (330)  796-4747. 

MIL-STD-2407  defines  the  VPF  standard.  VPF  data,  also  called  meta-data,  is  arranged  as 
directories,  tables  and  indices.  Essentially,  VPF  provides  a  model  that  describes  the  structure, 
organization,  and  relationships  of  the  information  in  a  terrain  database.  This  allows  fast,  direct 
access  of  the  database  with  no  need  for  translation.  The  NIMA  web  site  at  http://www.nima.mil 
has  more  information  about  VPF. 

SEDRIS  uses  the  concept  of  meta-data  and  extends  it  to  the  synthetic  environment.  SEDRIS 
calls  for  an  application  programmer  interface  (API)  to  access  a  terrain  database.  An  API  converts 
between  a  simulation’s  native  data  format  and  the  SEDRIS  model.  This  means  that  simulations  in 
a  synthetic  environment  can  truly  be  interoperable.  Instead  of  relying  on  custom  terrain 
databases,  each  simulation  in  the  synthetic  environment  could  simply  use  an  API  to  interact  with  a 
common  terrain  database.  Another  distinct  advantage  of  the  SEDRIS  model  is  that  it  is  easier  to 
communicate  with  other  simulations  by  describing  the  data  through  attributes  than  through  a  data 
storage  format.  For  more  information,  visit  the  SEDRIS  home  page  on  the  web  at 
http://www.sedris.net. 

6.0  Conclusion 

Terrain  database  development  is  labor  intensive  and  time-consuming.  However,  a  GIS  manager 
can  organize  a  tool  chest  of  software  that  makes  building  terrain  databases  easier.  Careful 
planning  and  ample  resource  allocation  ensures  that  the  terrain  database  is  completed  quickly. 
Once  the  terrain  database  has  been  developed,  have  all  the  simulations  in  the  synthetic 
environment  use  a  common  terrain  database.  By  eliminating  the  need  for  terrain  database 
conversion  between  simulations,  you  speed  up  data  access  and  minimize  the  map  errors  that 
inevitably  plague  distributed  simulations. 
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Appendix  A  -  Feature  Identification  Codes  supported  by  ARIES 


The  following  Feature  Identification  (FID)  codes  are  used  to  describe  the  predominant  nature  of 
all  features  (area,  linear,  and  point)  that  are  supported  by  ARIES. 

Feature  Identification  FID  Code 


Area  Features 


Quarry . 102 

Depot . 778 

Soil . 902 

Packed  Sand  &  Gravel . 906 

Sand  Dunes . 907 

Salt  Marsh . 908 

Smooth  Solid  Rock. . 910 

Rocky  Flat . 912 

Dry  Lake . 913 

Flood  Plain . 914 

Loose  Sand . 917 

Dry  Depression . 918 

Wadi . 919 

Salt  Flat . 934 

Fresh  Water  (General) . 940 

Non-Perennial  Stream  (Linear  Portrayal.) . 945 

Orchard/Hedgerow  (Background) . 951 

Irrigated  Field . 958 


Linear  Features 


Railroad . 206 

Dual  Highway  (with  Mediaa) . 250 

All  Weather  Hard  Surface  Highway . 251 

All  Weather  Loose  or  Light  Surface  Road . 252 

Fair  Weather  Loose  or  Light  Surface  Road . 253 

Cart  Track,  Trail . 254 

Road,  Approximate  Alignment,  Under  Construction, 

Existence  Reported . 255 

Pipeline  (Above  Ground). . 281 

Powerline  Pylon  (Type  “A”) . 541 

Powerline  Pylon  (Type  “H”) . 542 

Powerline  Pylon  (Type ‘T) . 543 

Powerline  Pylon  (Type  “Y”) . 544 

Runway  and  Taxi  way . 706 

Cleared  Way. . 916 

Wadi . 919 

Levee . 921 

Wall . 922 

Escarpment . 924 

Chain  Link  Fence . 927 

Fresh  Water . 940 

Non-Perennial  Stream  (Linear  Portrayal) . 945 

Canal/Channelized  Stream/Drainage  Ditch, 

(Subject  to  Ice,  Linear  Portrayal) . 947 

Revetment . 981 

Berm  . 982 
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Barbed  Wire  Fence . 983 

Concertina  Fence . 984 

Ditch . 985 

Trench . 986 


Point  Features 

Gas/Oil  Derrick . : . 103 

Offshore  Platform . 104 

Refinery . 120 

Power  Plant  (General) . 130 

Substation. . 138 

Light  Fabrication  Industry  (General) . 160 

Associated  Structure  (General  Industry.) . 180 

Building . 181 

Smokestack . 182 

Pumping  Station . 184 

Railroad  Station. . 222 

Bridges  (General) . 260 

Associated  Structure  (General  Transportation) . 290 

Commercial  Building  (General) . 301 

Grandstand . 324 

Multi  “Family  Dwelling  (General) . 401 

Single  Family  Dwelling  (General) . 420 

Agricultural  Building  (General) . 430 

Cemetery  Building . 451 

Communication  Tower . 501 

Miscellaneous  Tower . 530 

Power  Transmission  Tower  (General) . 540 

Governmental  (General) . 601 

Prison . 604 

School . 620 

Hospital . 630 

House  of  Religious  Worship  (General) . 650 

Associated  Structure  (General  Institutional) . 680 

Airport/Airbase  (General) . 701 

Ground  Support  Facility  (General) . 770 

Tank  (General) . 801 

Grain  Elevator . 822 

Water  Tower  (Building) . 824 

Warehouse . 861 

Date  Palm . 957 
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Appendix  B  -  Modified  DFAD  Feature  Attributes 


FID 

HEIGHT 

SMCCODE 

ORIENTATION 

DIRECT 

LENGTH 

WIDTH 

(m) 

(deg) 

(m) 

(m) 

251 

0 

14 

- 

3 

- 

4 

252 

0 

14 

- 

3 

- 

4 

253 

0 

5 

- 

3 

- 

4 

ISM 

0 

5 

- 

3 

- 

2 

2 

2 

- 

2 

- 

1 

706 

0 

9 

- 

3 

- 

25/12/8/4 

921 

3/1 

5 

- 

2 

- 

6/3 

922 

3 

3 

- 

2 

- 

2 

927 

3 

1 

- 

2 

- 

1 

940 

0 

6 

- 

2 

- 

6 

947 

0 

6 

- 

2 

- 

15 

984 

2 

1 

- 

2 

- 

1 

15 

2 

0-8 

- 

25 

25 

15 

2 

0-8 

- 

25 

25 

138 

5 

3 

0-8 

- 

15 

15 

180 

5 

3 

0-8 

- 

15 

15 

181 

5 

3 

0-8 

- 

15 

15 

182 

15 

3 

0-8 

- 

2 

0 

184 

5 

3 

0-8 

- 

15 

15 

222 

4 

2 

0-8 

- 

25 

25 

420 

3/2 

3 

0-8 

10/6 

4 

5 

3 

0-8 

- 

15 

15 

5 

3 

0-8 

- 

15 

15 

650 

5 

3 

0-8 

- 

15 

10 

701 

3 

3 

0-8 

- 

12/5 

12/5 

770 

2 

3 

0-8 

- 

10/2 

8/2 

■mol 

4 

1 

0-8 

- 

4 

0 

6 

1 

0-8 

- 

2 

0 

5 

3 

0-8 

- 

25 

25 

■SB 

4 

12 

0-8 

- 

2 

0 
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Appendix  C  -  Digital  Data  Base  ICD  for  ARIES 


DIGITAL  DATA  BASE 
Interface  Control  Document 

for  the 

ADVANCED  RADAR  IMAGING 
EMULATION  SYSTEM  (ARIES) 


Contract  Number:  F33615-95-C-1610 
CDRL  Sequence  Number:  A025 
30  September  1996 
Prepared  for: 

Air  Force  Wright  Laboratory 
Prepared  by: 

Lockheed  Martin  Tactical  Defense  Systems 
Post  Office  Box  85 
Litchfield  Park,  Arizona  85340-0085 


This  document  includes  proprietary  data  that  shall  not  be  disclosed  outside  Lockheed  Martin 
Tactical  Defense  Systems  and  its  designated  parties  and  shall  not  be  duplicated,  used,  or 
disclosed,  in  whole,  or  in  part,  for  any  purpose. 

This  restriction  does  not  limit  the  customer’s  right  to  use  information  contained  in  the  data  if  it  is 
obtained  from  another  source  without  restriction. 
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1.  SCOPE 

This  document  defines  the  contents  of  the  digital  terrain  and  feature  data  base  to  be  used  by  the 
ARIES  Synthetic  Aperture  Radar  (SAR)  Imagery  Simulation  being  developed  by  Lockheed- 
Martin  Tactical  Defense  Systems  (Lockheed-Martin)  for  incorporation  into  the  Radar  Processor 
Simulation  (RPS)  on  the  Joint  STARS  platform.  The  RPS  is  being  developed  for  the  integration 
of  Joint  STARS  into  the  Joint  Advanced  Distributed  Simulation  (JADS)  environment. 

1.1  Purpose 

The  purpose  of  this  document  is  to  define  the  characteristics  and  structure  of  the  digital  terrain 
elevation  and  feature  data  base  for  the  ARIES  simulation.  The  SAR  image  produced  by  the 
ARIES  simulation  represents  the  terrain  elevation  characteristics  and  the  specific  features  and 
their  locations  in  the  area  simulated.  This  data  base  will  be  structured  Defense  Mapping  Agency 
Digital  Terrain  Elevation  Data  and  Digital  Feature  Analysis  Data  products. 

1.2  Application 

Interface  requirements  set  forth  in  this  document  apply  during  the  development  and  testing  of  the 
ARIES  SAR  simulation  and  the  RPS. 

1.3  Definitions  and  Conventions 

The  following  conventions  were  used  to  describe  each  message  interface: 


2.  Applicable  Documents 

The  following  documents  are  applicable  to  the  extent  specified  herein: 

DMA  Digital  Terrain  Elevation  Data  (DTED)  Specification 
DMA  Digital  Feature  Analysis  Data  (DFAD)  Specification 

3.  Data  base  descriptions 

3.1  Data  Base  Coordinate  System 

The  data  base  shall  utilize  a  topocentric  coordinate  system.  The  coordinate  system  uses  a 
reference  plane  tangent  to  the  earth  at  the  Latitude  and  Longitude  specified  when  the  data  base  is 
constructed.  Lines  in  the  reference  plane  are  orthogonal  while  lines  of  Longitude  on  the  earth’s 
surface  curve  together  as  Latitude  increases. 

3.1.1  Units  of  Measure 

All  measurements  shall  be  in  meters. 

3.1.2  Horizontal  Reference 

The  axes  in  the  plane  shall  be  oriented  East-West  (X)  and  North-South  (Y).  Displacements  to  the 
North  and  East  from  the  topocentric  center  shall  be  positive.  Displacements  to  the  South  and 
West  shall  be  negative. 

3.1.3  Vertical  Reference 

Elevation  at  the  data  points  shall  be  referenced  to  the  topocentric  plane.  Points  with  elevations 
below  the  plane  shall  be  negative,  points  above  the  plane  shall  be  positive.  Points  which  fall  in  the 
plane  will  have  an  elevation  of  zero.  It  should  be  noted  that  points  of  constant  elevation 
referenced  to  Sea  Level  in  the  DTED  data  base  will  produce  elevations  that  vary  in  the 
topocentric  coordinate  system,  based  on  the  distance  from  the  JSTARS  topocentric  center. 
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3.1.4  Maximum  Coordinates 

The  maximum  displacement  along  either  horizontal  axis  from  the  topocentric  center  shall  be  ± 
256,000  meters.  The  maximum  displacement  along  the  vertical  axis  shall  be  -12,000  meters  to 
+9,000  meters. 

3.2  Digital  Terrain  Data  Base 

The  ARIES  digital  terrain  data  base  will  be  derived  from  DMA  DTED  Level  1  terrain  elevation 
data.  The  terrain  data  base  will  merge  data  from  a  number  of  standard  DMA  DTED  files  to 
produce  one  data  base  covering  the  entire  512  KM  x  512  KM  area. 

A  coordinate  conversion  will  be  required  to  map  the  DTED  data  referenced  to  a  Latitude- 
Longitude  coordinate  system  to  the  JSTARS  topocentric  coordinate  system. 

3.2.1  Terrain  Database  Structure 

The  terrain  database  shall  be  structured  in  the  same  manner  as  specified  in  the  DTED 
specification.  The  number  of  entries  at  any  given  latitude  will  be  a  constant  due  to  the 
orthogonality  of  the  coordinate  system  compared  to  the  Latitude-Longitude  coordinate  system. 

3.3  Digital  Feature  Data  Base 

The  ARIES  digital  feature  data  base  will  be  derived  from  DMA  DEAD  Level  1  feature  data.  The 
feature  data  base  will  merge  data  from  a  number  of  standard  DMA  DEAD  files  to  produce  one 
data  base  covering  the  entire  512  KM  x  512  KM  area. 

A  coordinate  conversion  will  be  required  to  map  the  DEAD  data  referenced  to  a  Latitude- 
Longitude  coordinate  system  to  the  JSTARS  topocentric  coordinate  system. 

Additional  features  may  be  added  to  the  data  base  at  the  direction  of  the  IADS  program  office. 
Eormat  for  these  features  shall  be  in  accordance  with  the  DMA  DEAD  specification.  New  types 
of  features  will  also  be  added  which  are  not  represented  by  current  feature  Identification 
Numbers  (EID). 

The  table  below  defines  the  new  features  (non-Df  AD)  to  be  added  to  the  Digital  feature 
Database. 

Non-DFAD  Feature  Definitions 


NAME 

FID 

TYPE 

Chain  Link  Fence 

935 

linear 

Barbed  Wire  Fence 

936 

linear 

Concertina  Fence 

937 

linear 

Anti-Tank  Ditches 

938 

linear 

Date  Palm  Trees 

957 

point 

Some  features  required  are  not  directly  supported  by  DEAD  features  but  can  be  indirectly 
supported  using  corresponding  existing  DEAD  f  ID’s.  The  table  below  identifies  these  features 
and  the  substitute  feature  to  be  used. 
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Substitute  DFAD  Feature  Definitions 


Non  DFAD  Feature 

Substitute  DFAD  Feature 

Rocky  flats 

Boulder  Field  (FID(#91 1),  Roeky,  Rough  Surface  (FID  #912) 

Packed  sand  and  gravel 

Sand/Desert  (FID  #906) 

Loose  sand 

Sand/Desert  (FID  #906) 

Dry  depressions  with  sandy  bottoms 

Sand/Desert  (FID  #906)  with  DTED 

Wadis 

Sand/Desert  (FID  #906)  with  DTED 

Escarpments 

Ground  Surface  (FID  #902),  Sand/Desert  (FID  #906),  Cliffs  (FID 
#924)  with  DTED 

Salt  marshes 

Ground  Surface  (FID  #902),  Marsh/Swamp  (FID  #908) 

Salt  flats 

Salt  Pans  (FID  #934) 

Flood  plains 

Ground  Surface  (FID  #902),  Mud/Tidal  Flats  (FID  #914) 

Date  palm  orchards 

Orchards  (FID  #951) 

Irrigated  fields 

Soil  (FID  #902),  Vegetation  (FID  #950) 

Oil  wells 

Gas/Oil  Derrick  (FID  #103) 

Revetments 

Levees/Embankments  (FID  #921),  Low  Embankments/Low  Levees 
(FID  #980) 

Below  ground  sand/dirt  trenches 

Ground  Surface  (FID  #902),  Sand/Desert  (FID  #906)  with  DTED 

Sand/dirt  ditches 

Ground  Surface  (FID  #902),  Sand/Desert  (FID  #906)  with  DTED 

Transmission  towers  -  4  sided 
pyramidal 

Communications  Towers  (FID  #501),  Radio/Television  Towers  (FID 
#’s  51 1,  512),  Power  Transmission  Towers  (FID#504) 

Electrical  power  lines 

Powerline  Pylons  (FID  #’s  541-544) 

Dirt  and  concrete  dikes  and  levees 

Conduits  (FID  #280),  Levees/Embankments  (FID  #921),  Low 
Embankments/Low  Levees  (FID  #980) 

Dirt  and  concrete  walls  and  berms 

Walls  (FID  #922) 

Pipelines  within  trenches 

Pipelines  (Above  Ground)  (FID  #281) 

3.3.1  Feature  Data  Base  Structure 

The  feature  database  shall  be  structured  in  the  same  manner  as  specified  in  the  DMA  DFAD 
specification.  Changes  to  this  structure  will  occur  in  the  maximum  number  of  features  per  Hgt?) 
set  and  the  limitations  on  the  maximum  geographic  area  covered  by  the  data  base. 

4.3  Digital  Contour  Terrain  Data  Base 

The  ARIES  digital  contour  terrain  data  base  will  be  derived  from  DMA  DTED  Level  1  terrain 
elevation  data.  The  contour  terrain  data  base  will  merge  data  from  a  number  of  standard  DMA 
DTED  files  to  produce  one  data  base  covering  the  entire  512  KM  x  512  KM  area. 

A  coordinate  conversion  will  be  required  to  map  the  DTED  data  referenced  to  a  Latitude- 
Longitude  coordinate  system  to  the  JSTARS  topocentric  coordinate  system. 
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4.3.1  CONTOUR  TERRAIN  Data  Base  Structure 

The  contour  terrain  database  shall  be  structured  in  contour  vectors.  Each  contour  vector  shall  be 
separated  in  elevation  by  10  meters.  The  contour  terrain  database  file  will  a  binary  data  file  as 
created  by  PV-WAVE®. 
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Appendix  D  -  ARC/INFO®  Macro  Language 
Scripts 


/*  startup. ami 

/*  Unix  version  of  startup. ami 

/♦ 

/*  Run  this  ami  to  set  up  the  Arc/Info  environment  and  initialize 
/*  convenient  variables  for  each  workspace. 

/* 

/*  Created  1/2/97  by  Capt  Rodney  Houser 

/■* 

&tcrm  9999 
display  9999  3 
coordinate  mouse 

&sctvar  .etc  = /diskl/uscrs/ai/ctc/ 

&sctvar  .arics  =  /diskl/uscrs/ai/ctc/arics/ 

&sctvar  .janus  =  /disk  1  /uscrs/ai/ctc/sccnario/janu.s/ 
precision  double 
&  return 


/*  join_codcs.aml 

/*  This  ami  joins  the  pcodc,  acodc,  and  xcodc  files  for  one  coverage.  Copy 
/*  and  run  this  ami  from  the  workspace  that  contains  the  coverage. 

/* 

/*  Created  08/07/97  by  Capt  Rodney  Houser 
/* 

relate  restore  dsOlp.relatc 
tables 

additem  dsOIp.aat  fucodc  4  5  B  #  dsOIp-id 

scl  dsOIp.aat 

calc  facodc  =  dsOlp-id 

q 

joinitem  dsOIp.aat  dsOIp.pcodc  dsOIp.aat  facodc  facodc 
rcgionclass  d.s01p  reg  land  dsOlp-id  facodc 
clean  reg  reg  #  H  poly 

joinitem  rcg.patland  rcg.pcode  rcg.patland  facodc  facodc 

crcatclabcis  reg 

relate  restore  dsOl  1. relate 

tables 

additem  dsOll.aat  facodc  4  5  B  #  dsOll-id 

scl  dsOll.aat 

calc  facodc  =  ds01I-id 

q 

joinitem  dsOll.aat  dsOll.acodc  dsOll.aat  facodc  facodc 
tables 

additem  dsOl  l.pat  facodc  4  5  B  #  dsOl  1-id 

scl  d.s01  l.pat 

calc  facodc  =  dsO  11-id 

q 

joinitem  dsOll.pat  dsOll.xcodc  dsOlI.pat  facodc  facodc 
&  return 


/*  crcatc_rcgions.aml 

/*  This  ami  creates  a  new  region  subclass  LAND  for  DFADARC  coverages,  and 
/*  joins  attribute  information  to  that  subclass.  Replace  XXXX 
/*  with  the  workspace  variable  name. 

I* 

/*  Created  08/07/97  by  Capt  Rodney  Houser 

/• 

relate  restore  dsOlpXX XX. relate 
tables 

select  dsOlpXXXX.aat 

additem  dsOlpXXXX.aat  facodc  4  5  B  #  dsOlpXXXX-id 
calc  facodc  =  ds01pXXXX-id 

q 

joinitem  dsOlpXXXX.aat  dsOlpXXXX.pcodc  dsOlpXXXX.aat  facodc  facodc 
rcgionclass  d.s01pXXXX  rcgXXXX  land  dsOlpXXXX-id  facodc 
clean  rcgXXXX  rcgXXXX  H  U  poly 

joinitem  rcgXXXX.patland  rcgXXXX.pcodc  rcgXXXX.patland  facodc  facodc 

crcatclabcis  rcgXXXX 

&rctum 


additem  dsOIp.pcodc  pin  4  5  B  #  facodc 

/* 

scl  dsOIp.pcodc 
calc  pin  =  Srccno 
scl  dsOll.acodc 
calc  pin  =  Srccno 
scl  dsOll.xcodc 
calc  pin  “  Sreeno 
scl 
I* 

dir  *codc 

&pausc  Identify  number  of  records 

&s  pcodc  [response  'Enter  number  of  pcodc  records'] 

&s  acodc  [response  'Enter  number  of  acodc  records'] 

scl  dsOll.acodc 

calc  pin  =  pin  +  %pcodc% 

scl  dsOll.xcodc 

calc  pin  =  pin  +  %pcodc%  +  %acodc% 

q 

/* 

pullitcms  dsOIp.pcodc  dsOIp.ppin 

~facodc 

~pin 

~cnd 

pullitcms  dsOll.acodc  dsOll.apin 

-facodc 

-pin 

-end 

pullitcms  dsOll.xcodc  dsOll.xpin 

-facodc 

-pin 

-end 

/* 

joinitem  dsOIp.aat  dsOIp.ppin  dsOIp.aat  facodc  facodc 
joinitem  dsOll.aat  dsOll.apin  dsOll.aat  facodc  facodc 
joinitem  dsOll.pat  dsOll.xpin  dsOll.pat  facodc  facodc 

/* 

tables 

scl  dsOIp.aat 
calc  ds01p-id  =  pin 
calc  facodc  =  pin 
scl  dsOll.aat 
calc  ds01l-id  =  pin 
calc  facodc  =  pin 
scl  dsOll.pat 
calc  ds01l-id  =  pin 
calc  facodc  =  pin 
scl 
/* 

dropitem  dsOl  l.pat  pin 
dropitem  dsOl  l.aat  pin 
dropitem  dsOIp.aat  pin 

I* 

scl  dsOIp.pcodc 
calc  facodc  =  pin 
scl 

dropitem  dsOl  p.pcodc  pin 
scl  dsOll.acodc 
calc  facodc  =  pin 
scl 

dropitem  dsO  1 1. acodc  pin 
scl  dsOll.xcodc 
calc  facodc  =  pin 
scl 

dropitem  dsOl  1. xcodc  pin 
/* 

kill  dsOIp.ppin 
kill  dsOll.apin 
kill  dsOll.xpin 

q 

/* 

idedit  dsOlp  line 
idedit  dsOll  line 
idedit  dsOll  point 
&rctum 


/*  modify _codcs.aml 

/*  This  ami  organizes  pcodc,  acodc,  and  xcodc  files  so  that  features 
/*  arc  incremented  by  1  in  ascending  numerical  order.  Copy  and  run 
/*  this  ami  from  the  workspace  that  contains  the  final  dsOlp  and 
/*  dsOll  coverages. 

/* 

/*  Created  08/07/97  by  Capt  Rodney  Houser 

/* 

tables 

additem  dsOll.acodc  pin  4  5  B  #  facodc 
additem  dsOl  l.xcodc  pin  4  5  B  #  facodc 
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Appendix  E  -  C  Code 


int  filctypcjd ; 


i^includc  <sys/filc.h> 

int  iopcn  (ch,  nch) 
char  *ch; 
int  *nch; 

{ 

return  (open  (ch,  0_CREAT  |  O.WRONLY,  0777) ); 

} 

int  iopcnr (ch,  nch) 
char  *ch; 
int  *nch; 

{ 

return  (open  (ch,  0_RD0NLY)  ); 

} 

int  ircad  (fd,  buf,  nbytes) 
int  *fd; 
char  *buf; 
int  *nbytcs: 

{ 

int  num; 

num  =  read  (*fd,  buf,  *nbytcs); 
if  (num  <  0) 

{ 

penror  ('^nbinio  read"); 

} 

return  (num); 

} 

int  iwritc  (fd,  buf,  nbytes) 
int  *fd; 
char  *buf; 
int  *nbytcs: 

{ 

int  num; 

num  =  write  (*fd,  buf,  *nbytcs); 
if  (num  <  0) 

{ 

perror  ('^nbinio  write"); 

} 

return  (num); 

} 


declare  functions  ♦♦♦♦****************************'*"'‘**^***/ 

void  writc„fcaturc„sct() ; 

int  gcMokcns  (char  ^string,  char  *scparator,  char  *tokcn[], 
int  max_tokcn_chars,  int  max„num_tokcns) ; 


global  variables  *****************************************/ 

int  num^fcaturcs  =  0  ; 
int  fcaturc_num  ; 
int  npts ; 

double  fcaturc_casting  [MAX_NLfM_VERTICES_PER_FEATUREI ; 
double  fcaturc_northing[MAX_NLJM_VERTlCES_PER_FEATURE] ; 

FILE  *fout; 
char  filctypc  [256] ; 


MAIN  PROGRAM  ********************************************* 

void  main  Q 
{ 

FILE  ♦fin ; 
int  i ; 

int  good_filc_typc ; 

char  filename  [256] ; 

char  outputfilc  [256] ; 

char  filctypes  [NUM_FILETYPES][25] ; 

int  numjokens ; 
char  linc[256] ; 

char  scparator[NUM_SEPARATOR_CHARS]  =  ’’  \t"  ;  /*  SPACES  OR  TABS  */ 
char  ♦tokcn[MAX_NUM_TOK:ENS] ; 


//  ALLOCATE  STORAGE  FOR  THE  TOKENS 
for  (i=0;  i<MAX„NUM_TOKENS;  i++)  { 
if  ( (tokcn[i]  =  (char  *)  malloc  (MAX_TOKEN„CHARS))  =  NULL)  { 
fprintf  (stderr,  "Couldn't  allocate  token  storagc\n")  ; 

} 

} 


/*  DEFINE  THE  ACCEPTABLE  FILE  TYPES  V 

strepy  (filctypcs[0),  "Dual  Highway") ; 

strepy  (filctypcs[l],  "Highway") ; 

strepy  (filctypcs[2],  "Road") ; 

strepy  (filctypcs[3],  "River") ; 

strepy  (filctypcs[4],  "Runway")  ; 

strepy  (filctypcs[5],  "Levee") ; 

strepy  (filetypcs[6],  "Wall") ; 

strepy  (filctypes  [7],  "Fence") ; 

strepy  (filctypcs[8],  "Lake") ; 


/*  PROMPT  USER  FOR  INPUT  FILE  */ 
printf("\n%s\n%s",  PROGJl,  PROG_I2) ; 
printf  ("\n\nEntcr  file  to  convert:  ") ; 
scanf  ("%s",  filename) ; 
printf  ("File  Selected:  %s",  filename) ; 


/*  CREATE  THE  OUTPUT  FILENAME  V 

strepy  (outputfilc.  filename) ; 

strepy  (strrehr  (outputfilc,  'O+l,  "env")  ; 


/*  OPEN  THE  OUTPUT  FILE  V 
if  ( (fout  =  fopen  (outputfilc,  "w"))  ==  NULL)  ( 
fprintf  (stderr,  '^nERROR  ..  could'nt  open  ouput  file:  %s\n", 
outputfilc) ; 

exit  (0) ; 

} 


/♦  PROMPT  USER  FOR  THE  TYPE  OF  FILE  ♦/ 
printf  ("\n\nID  Source  file  type") ; 

for  (i=0;  i<NUM„FILETYPES;  i++)  printf  ("\n%d  %s",  i,  filctypcs[i]) ; 

good_filc_typc  =  FALSE ; 
while  (!good_filc_typc)  { 

printf  ('^n\nEntcr  ID  of  Source  file  type  (99  to  exit):  ")  ; 
scanf  ("%d",  &filctypc_id) ; 

if  (filctypc_id  >  -1  &&  filctypc_id  <=  NUM_FILETYPES)  { 
good_filc_typc  =  TRUE ; 

} 

if  (filctypcjd  ==  99)  { 

printf  OnXnExiting  ..  thanks,  its  been  fun?!!\n")  ; 
exit  (0) ; 

} 

if  (!good_filc_typc)  { 

printf  ('^nYou  must  enter  the  ID  of  one  of  the  Source  File  types") ; 
printf  ("  or  99  to  exit !!!")  ; 

} 

) 

printf  ("ID  of  Source  file  entered:  %d  Filctypc:  %s\n". 


/*  arc2janus.c  */ 

ffdcftnc  PROGJl  "This  program  converts  an  arcinfo  ungcncratc  ascii  file" 
/(define  PROGJ2  "  (in  utm)  into  janus  compatible  terrain  data  input  file" 

//include  <stdio.h> 

//include  <string.h> 

//include  <stdlib.h> 


//define  NUM_FILETYPES  9 

//define  TRUE  1 

//define  FALSE  0 

^define  MAX_NUM_TOKENS  5 

//define  MAX_TOKEN_CHARS  20 

//define  NUM„SEPARATOR_CHARS  4 

//define  MAX_NUM_FEATURES  10000/*  janus  limit  is  10000  */ 

//define  MAX„NUM_VERTICES  100000  /*  janus  limit  is  100000  */ 

//define  M AX_NUM_VERTICES_PER_FEATURE  1 000  /*  janus  limit  is  1000  */ 
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filctypcjd,  filctypcs[filctypc_id]) ; 


I*  OPEN  THE  INPUT  FILE  */ 
if  ( (fin  =  fopcn  (filename,  "r"))  —  NULL)  { 
fprintf  (stderr,  "\nUnablc  to  open  %s\n",  filename) ; 

) 


/*  PERFORM  THE  CONVERSIONS  */ 
while  (fgets  (line,  sizcof  (line),  fin)  !=  NULL)  { 

*strchr  (line,  W)  =  ^0' ;  II  REMOVE  NEWLINE 

numjokens  =  gct_tokens  (line,  separator,  token. 

MAX_TOKEN_CHARS,  MAX„NUM_TOKENS) ; 

if  (num_tokens  ==  0)  { 

fprintf  (.stderr,  "\nERROR  ..  got  zero  tokcns\n") ; 
break ; 

} 


//  IF  ONLY  1  TOKEN,  IT'S  THE  END  OF  A  FEATURE  SET  OR  BEGINNING 
if  (num_tokcns  ==  1)  { 
if  (stremp  (tokenfO],  "END")  =  0)  { 
if  (npLs  ==  0)  { 

printf  ("\n\nTHArS  ALL  FOLKS  ..  new  file;  %s\n",  outputfilc) ; 
cxit(O) ; 

}  else  { 

writc_fcaturc_set  () ; 
npts  =  0 ; 

}  else  { 

fcaturc_num  =  atoi  (lokenfO]) ; 
if  (featurc_num  ==  0)  { 

fprintf  (.stderr,  "\nERROR  ..  got  feature  number  of  0")  ; 
fprintf  (stderr,  "\n  could  be  bad  data\n") ; 

exit  (0) ; 

} 

} 

} 

//  IF  2  TOKENS,  GET  THE  EASTING,  NORTHING  FEATURE  DATA 
if  (nuni .tokens  “  2)  { 
fcaturc.casting  [npts]  =  atof  (tokcn[0]) ; 
fcaturc_northing(npts]  =  atof  (tokcn[l]) ; 
npt.s++ ; 

if(npts>  MAX_NUM_VERTICES_PER_FEATURE)  { 
fprintf  (stderr,  "\nERROR  ..  exceeded  max  number  vertices'^") ; 
exit  (0) ; 

} 

) 

) 

}  /*  END  OF  MAIN  */ 


/♦*♦*♦*♦♦*  WRITE_FEATURE_SET***************^***********’^************/ 
void  write_fcaturc_sct  Q 
{ 

int  i ; 

char  hcadcr.str  [MAX_NUM_FEATURES][256] ; 

strepy  (hcadcr_str[0],  "6\\1\\1\\  16.\\Dual  Highway") ; 
strepy  (hcadcr_str[l],  "6\\2\\2\\  8,\\Highway") ; 
strepy  (headcr_str[2],  "6\\3\\3\\  S.WRoad") ; 

.strepy  (hcadcr_str[3],  "2\\1\\1\\  lOO.WRivcr")  ; 
strepy  (hcadcr_str[4],  "3\\1\\1\\  B.WRunway") ; 
strepy  (hcadcr_str[5),  "3\\2\\2\\  12A\Lcvcc")  ; 
strepy  (hcadcr.str[6],  "3\\3\\3\\  l.WWall") ; 
strepy  (hcadcr_str[7],  "4\\1\\1\\  O.WFcncc") ; 
strepy  (hcadcr_str[8],  "5\\l\\l\\0A\Lakc") ; 


GET  TOKENS  ****’<'***************’*‘**************************/ 
int  gct.tokcns  (char  *string, 
char  *  separator, 
char  '*tokcn[], 
int  max_tokcn_chars, 
int  max_num_tokcns) 

{ 

/*  obtains  tokens  in  string  *! 

I*  returns  number  of  tokens  found,  0  if  unsuccessful  *1 

I*  *1 

I*  the  following  must  be  defined  prior  to  calling  */ 

I*  token:  pointer  storage  must  be  allocated  *1 
I*  separator:  defines  the  chars  that  separate  tokens*/ 

/*  max_tokcn_chars;  max  num  of  chars  allowed  in  a  token  */ 

/*  max.num.tokcns:  max  num  of  tokens  */ 

y  +  +  **  +  «***  +  **1>*  +  *  +  *  +  i»K>*+i|>KiiMHC*  +  **  +  ***  +  1>*>ti*>l"***t>*f*  +  +  ***'t“l‘**/ 


char  *.strptr=" 

int  Icn,  num.tokcns  =  0  ; 


II  MAKE  SURE  THE  INPUT  IS  NOT  A  NULL  STRING 
if  (  ! stremp  (string, ""))  return  (0) ; 


n  PULL  OFF  THE  TOKENS 
while  (strptr  !=  NULL)  { 
if  (numjokens  ==  0)  { //  GET  FIRST  TOKEN 
strptr  =  strtok  (string,  separator) ; 
if  ( (ten  =  strlcn(strptr))  >  max_tokcn_chars)  { 
fprintf  (stderr,  "\nERROR..gct_tokcns") ; 
fprintf  (stderr.  'An  first  token  had  %d  chars",  Icn) ; 
fprintf  (stderr,  "\n  max  token  chars: ") ; 
fprintf  (stderr,  "%d",  max_tokcn_chars); 
fprintf  (stderr,  'An") ; 
return  (0) ; 

}  else  { 

strepy  (tokcn[num_tokcns++],  strptr) ; 

} 

)clsc{  //OTHER  TOKENS 

strptr  =  strtok  (NULL,  separator) ; 
if  (strptr  !=  NULL)  { 

If  (num_tokcns  >  max_num_tokcns)  { 
fprintf  (stderr,  'AnERROR..gct_tokcns") ; 
fprintf  (stderr,  'An  exceeded  max  num  tokens: ") ; 
fyrintf  (stderr,  "%d",  max_num_tokcns): 
fprintf  (stderr,  "\n") ; 
return  (0) ; 

}  else  { 

if  ( (Icn  "  strlcn(strptr))  >  max_token_char.s)  { 
fprintf  (stderr,  "\nERROR  ..  gct.tokens")  ; 
fprintf  (stderr,  "\n  subs  token  had  %d  ehars",  Icn); 

fprintf  (stderr,  "\n  max  token  chars: ") ; 
fprintf  (stderr,  "%d",  max_tokcn_chars)  ; 
fjjrintf  (stderr,  'An") ; 

}  else  { 

strepy  (tokcn[num_tokcns++],  strptr) ; 

) 

} 

} 

) 

} 

return  (num_tokcns) ; 

}  /*  END  OF  GET.TOKENS  */ 


printf  ("\nwriting  fcaturc.num;  %d  with  %d  vertices",  feature.num,  npts) ; 

fprintf  (fout,  "%d  %s  %I0d\n", 

fcaturc.num,  headcr_str[filctype_id],  npts) ; 

for  (i=0;  i<npts;  i++) 
fprintf  (fout,  "%  1 0.2f  %  1 0.2 An", 

fcaturc_easting[i],  fcaturc_northing[i]) ; 

}  /*  END  OF  WR1TE_FEATURE_SET  */ 
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Appendix  F  -  FORTRAN  Code 


PROGRAM  DFAD 
C 

CSUBR0UTINES:TAPEI0J,TAPEI0_C,C32T35,B1N2DEC_D,WGSUTM, 

C  BlN2DEC„A,INOUT,PTGEN.PLOT:PLOT/LIB,UTMWGS 
C  LINKED  BY  :link.dfad.9.com 

C  R  INTERGER  MATRIX  HOLDING  THE  FIC  NUMBER  OF  THE  ROADS 
C  V  BYTE  MATRIC  HOLDING  THE  SURFACE  FEATURE  HEIGHTS 

COMMON/A/BUFF36(I2.8),BUFFBITS(288) 


INTEGER*4  SBUFF(9),ifiIc 

CHARACTER*80  OUTFILE,infilc,filcnanic 
byte  null 

character*!  anull,<u;pacc 
BYTE  BUFF36,BUFFBITS 
cquivalcncc(null,anull) 


C************  READ  INPUT  PARAMETERS  ♦****♦♦**♦*•***♦*3 

aspacc  = ' ' 
null  =  0 

writc(6,"CEntcr  in  DFAD  file  to  process  (no  ext)  ==>',$)") 

rcacl(5,'(a40)')filcfi3mc 

idx  =  indcx(filcnamc, aspacc)-! 

infilc  =  fiIcnamc(l:idx)/AlvcV/anuII 

outfilc  =  filcnamc(l:idx)//’.36bits’ 

print  *,infilc 
print  *, outfilc 

ifilc=  iopcnifinfilc) 

OPEN  (UNIT-3, FILE=OUTFILE, 

*  STATUS=’UNKNOWN’, 

*  FORM='FORMATTED') 

lRCNT-1 
ic  =  0 
itot  =  0 

100  CONTINUE 
ICOUNT=0 


num  =  ircad  (ifilc,SBUFF,36) 
irent  =  irent  +  1 
itot  -  itot  +  num 
if  (mod(ircnt,100).cq.0)thcn 
print  *,'Rcc,bytcs  read, total  bytes  ',ircnt,num,itot 
endif 

if(num.cq.0)go  to  999 
c  if(ic.gt.l0000)go  to  999 

CALL  C32T36(SBUFF) 

do  ilp-  1,8 
ic  =  ic  1 

writc(3,l  1 1)(buffbits(ij),ij=C36*(ilp-l)+l),(36*ilp)).ic 
Ill  fonnat(36il,5x,il5) 
enddo 

do  ij=l  ,9 
sbufRij)  =  0 
enddo 
do  ij=  1,288 
bunbits(ij)  =  0 
enddo 

if(l.cq.l)go  to  100 


999  CONTINUE 
CLOSE  (UNIT=3) 
STOP 
END 


******•¥********  *************** 

PROGRAM  DFAD 


C 

C  SUBR0UT1NES:B1N2DEC_D, 

CLINKED  BY  ilink.arics 

Q  ******************************************************** 

c 

C  This  program  will  read  DFAD  format  off  of  tape  and  processes 
C  data  into  a  .tvg  format 
C 

COMMON/A/BUFF36(12,8),BUFFBITS(288) 

COMMON/D/ADATA(21600*30) 


INTEGER*4  STATUS 

INTEGER*4IFEAID(1000,3),rSMCAR(14,3) 

1NTEGER*4  XINTNUM,ICOUNT,WAGWACN,WAGWACC,WAGCELL 
INTEGER*4  poricnt,plcngth,pwidth 
INTEGER*4  Idirccflwidth 
INTEGER*4  astruct,atrcc,aroof 

CHARACTER*80  OUTFILE,INFILE,filcnamc 
CHARACTER*!  DSI(648),ACC(2700),SPACE,anull 

BYTE  BlTS36(36),nulI 

BYTE  BUFF36.BUFFB1TS,ADATA 

DIMENSION  X(8000),Y(8000) 

EQUIVALENCE  (STAT.STATUS) 
equivalence  (null.anull) 

DATA  SPACE/’  ’/ 


Q************  INPUT  PARAMETERS  *****************s 

writc(6,"('Entcr  in  file  to  process  =>  ',$)") 

read(5,'(a40)')filcnamc 

idx  =  indcx(filcnamc,spacc)-l 

infile  =  filename(l  :idx)//'.36bits’ 

outfilc  =  filcnamc(l;idx)//'.dfad' 


Q  ****************  input  and  output  FILES  ♦♦*♦++*♦**♦**  + 

OPEN  (UNIT=2,FILE=INFILE, 

*  STATUS='UNKNOWN', 

*  FORM-’FORMATTED') 

OPEN  (UN1T=9,FILE=0UTFILE, 

*  STATUS-unknown', 

*  FORM='FORMATTED') 


doi=l,3 
doj-1,1000 
ifcaid(j,i)  =  0 
enddo 
enddo 


100  CONTINUE 
IRCNT  =1 
(COUNT  =  0 
NPTS  =0 
1REJECT  =  0 
IPHTMAX  =  0 
IHEAD  =0 
IDX  =  1 
ZONE  =38 

C - 

C-  Process  Manuscript  header.  Data  Set  Identification  Record,  Accuracy  Record 
C - 

c - 

c-  Manuscript  Data  Set  Header  consists  of  *6*  36-bit  words 

c - 

ice  =  1 
do  i=l,6 

rcad(2,500)bits36 
500  format(36iI) 
do  j=  1,36 

bufTbits(icc)  =  bits36(j) 
ice  =  lee  +  1 
enddo 
enddo 

DO  120  J=  1,288 
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ADATA(J)=BUFFBITS(J) 
120  CONTINUE 


CALL  BIN2DEC_D(1 ,6,XINTNUM) 

PRINT  *,XINTNUM 

C  IF  XINTNUM  =  63  — ->  END  OF  ALL  MANUSCRIPTS 
1F(XINTNUM.EQ.63)G0  TO  900 

!MAN=1MAN+1 
WRITE(6,I500)IMAN 
1500  F0RMAT(1X;  MANUSCRIPT  #M5) 

CALL  BIN2DEC_D(10.3,XINTNUM) 

ILEV=XINTNUM 

CALL  BIN2DEC_D(13, 14, XINTNUM) 
WAGWACN-XINTNUM 

CALL  BIN2DEC„D(27,5,XINTNUM) 
WAGWACC-XINTNUM 

CALL  BIN2DEC_D(32,5,XINTNUM) 
WAGCELL=X1NTNUM 

CALL  B1N2DEC„D(37,36,XINTNUM) 
ILAT10=XINTNUM 

CALL  BIN2DEC_D(73,36,X1NTNUM) 
ILON10=XINTNUM 

ALON=ILONl  0/1 073600. 

ALAT=]LAT1 0/1 073600. 
c  IF(ZONE.GT.30)ALON=-ALON 
c  WRITE(9,1605)ALAT,ALON 
1605  F0RMAT(1X;LAT,L0N’,2F12.4) 

CALL  B1N2DEC_D(109,18,XINTNUM) 
ILATMX=XINTNUM+ILAT10 

CALL  BIN2DEC_D(127.18, XINTNUM) 
1L0NMX=X1NTNUM+IL0N1 0 

ALONMX=ILONMX/l  073600. 

ALATMX=ILATMX/ 1 073600. 
c  IF  (ZONE.GT.30)ALONMX=-ALONMX 


WRITE(6.1100)IMAN.1LEV,WAGWACN,WAGWACC.WAGCELL,ALON,ALAT,ALONM 

X. 

♦ALATMX 


C - 

c-  DSI  consists  of  648  bytes  =  *144*  36-bit  words 


print  *,'Rcading  DSP 
DO  1=1,144 
READ(2.500)BITS36 
ENDDO 

c  WR1TE(9,1200)DSI(4),(DS1(I).1=7,33),(DSI(I).1=60,64). 
c  *(DSI(I),I=65,79),(DS1(I), 1=88,98), (DSia),I=127,141), 
c  *(DS1(I).1=145.149),(DSI(I),I=160,163) 

c-  ACC  consists  of  2700  bytes  =  *600*  36-bit  words 

print  *, 'Reading  ACC 
DO  1=1,600 
READ(2,500)BIT36 
ENDDO 

c  WRITE(9,1 300)(ACC(I).I=4.7),(ACC(I),1=1 2,15),(ACC(1),I=20,23), 
c  *  (ACC(I),I=56,57) 


1TOPLIM=21600 

C - 

C-  PROCESSING  FEATURES 
C - 

c  WRITE{9,2900) 

2900  FORMATC  ****  FEATURES  BEING  PROCESSED') 


i36  =  0 

300  CONTINUE 

C - 

C-  Read  Feature  data  header  -  *2*  36-BIT  WORDS 


do  i=l,288 
bufibits(i)  =  0 
enddo 

idx  =  l 
icc=  I 
do  i=l,2 

c - Checking  on  checksum  word- — 

ifl[i36.cq.600)thcn 
read(2,500)bits36 
c  writc(9,500)bits36 

rcad(2,500)bits36 
c  writc(9,500)bits36 

i36  =  0 

c  writc(9,*)'*****chccksum*****' 

endif 


READ(2,500)bits36 
i36  =  i36-)-l 

c  writc(9,*)i36,i36 
doj=l,36 

buflbits(icc)  =  bits36(j) 
ice  =  ice  +  1 
enddo 
enddo 

DOJ=l,288 

ADATA(J)=BUFFBITS(J) 

enddo 

c  writc(6, 500)  (buflbits(ij),ij=  1,72) 

IHEAD=IHEAD-H 

CALL  BIN2DEC_D{IDX,I4,XINTNUM) 

1FAC=XINTNUM 

CALL  BIN2DEC_D(IDX-H  4,2, XINTNUM) 

IFEATP=XINTNUM 

Q  ************H.***4.**+*4.*^  jEs-p  poR  end  of  MANUSCRIPT  ***+*■*■♦** 
IFaFEATP.EQ.3)GO  TO  800 
ITEST=IFEATP+1 

c  CALL  BIN2DEC_D(IDX+16+1.10-1,XINTNUM) 

CALL  BIN2DEC_DaDX-»-16,10,XINTNUM) 

IPHT=XINTNUM*2 

IF  (1PHT.GT.IPHTMAX)IPHTMAX=1PHT 

CALL  BIN2DEC_DaDX+26,lO,XINTNUM) 

1FICN=XINTNUM 

CALL  BlN2DEC_D(IDX+36,5, XINTNUM) 

ISMC=XINTNUM 

c - Point  Feature  Specifics - 

if(ifcatp.eq.O)thcn 

call  bin2dcc_d(idx+36+5,6,xintnum) 
porient  =  xintnum 

call  bin2dcc_d(idx+36+l  l,7,xmtnum) 
picngth  =  xintnum  *  2 

if  (ificn.ge.230.and.ificn.Ic.239)plcngth  =  xintnum  *  20 
call  bin2dcc_d(idx+36+l  8,7,xintnuni) 
pwidth  =  xintnum  *  2 

if  (ificn.gc.230.and.ificn.lc.239)pwidth  =  xintnum  *  20 
CALL  BlN2DEC_DaDX+61 ,1 1, XINTNUM) 

N=XINTNUM 

WR1TE(6,1400)IFAC,IFEATP.IPHT,IFICN,ISMC,N, 

*  poricnt,plcngth,pwidth 
endif 

if(ifcatp.cq.l)then 

call  bin2dec_d(idx-i-36-f-5,2,xintnum) 

Idircct  =  xintnum 

call  bin2dcc_d(idx+36+7,7, xintnum) 

Iwidth  =  xintnum  *  2 

cal  1  bin2dcc_d(idx-f  36+ 1 4, 1 4  ,xintnum) 

Iblank  =  xintnum 

CALL  BlN2DEC_DaDX+59,13,XINTNUM) 

N=XINTNUM 

WRITE(6,1401)IFAC,1FEATP,IPHT,!FICN,ISMC,N, 

*  ldircct,lwidth 
endif 

if{ifeatp.cq.2)then 

call  bin2dcc_d(idx+36+5,4,xintnum) 
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astruct  =  xintnum 

call  b  i  n2dcc_d(idx+3  6+9,4  ,x  intnum) 

atrcc  =  xintnum  *  10 

call  bin2dcc_d(idx+36+l  3.4, xintnum) 

aroof  =  xintnum  *  10 

call  bin2dcc_d(idx+36+1 7,6,xintnum) 

ablank  ~  xintnum 

CALL  BIN2DEC_D(IDX+59.13.XINTNUM) 
N=XINTNUM 

WRITE(6,1402)1FAC,1FEATP.1PHT.IFICN,ISMC,N, 
*  astruct.atrcc, aroof 

cndif 


IF  (IFICN.GT.O  .AND.  IFICN.LE.IOOO  .AND. 

*  ITEST.GT.O  .AND.  ITEST.LE.3) 

*  IFEAID(1F1CN,ITEST)=IFEAID(IFICN,ITEST)+I 
IF  (ISMC.GT.O  .AND.  ISMC.LE.I4  .AND. 

*  ITEST.GT.O  .AND.  ITEST.LE.3) 

*  ISMCARaSMC.lTEST)=ISMCARaSMC,ITEST)+l 
IF  (ISMC.LE.O  .OR.  ISMC.GT.14)]MISS=IMISS+1 


350  CONTINUE 

C - 

c-  Record  of  data  points  n#^  of  36-bit  words 

c - 

idx  =  l 
icc=  1 

c  print  *, 'Reading  ’,n,*  36-bit  words' 
do  i=l,n 

c - Checking  on  checksum  word - 

ifl[i36.cq.600)thcn 
rcad(2.500)bits36 
c  writc(9,500)bits36 

rcad(2.500)bits36 
c  writc(9,500)bits36 

i36  =  0 

c  writc(9,*)'*****Chccksum*****' 

cndif 

rcad(2,500)biis36 
i36  =  i36  +  1 
c  writc(9,*)i36 
doj=I,36 

adata(icc)  =  bits36(j) 
ice  =  ice  +  1 
enddo 
enddo 

NPTS=0 

DOJ==l,N 

NPTS=NPTS+1 

CALL  BIN2DEC_D(1DX+1 ,18-1, XINTNUM) 
c  Y(NPTS)=XlNTNUM/36000. 

Y(NPTS)=XINTNUM/36000.  +  ALAT 
c  if(y(npts).cq.alat.and.npts.gt.5)y(npts)  =  alatmx 

CALL  BIN2DEC_D(IDX+1 8+1, 18-1, XINTNUM) 
c  X(NPTS)=XINTNUM/36000. 

X(NPTS)=XlNTNUM/36000.  +  ALON 
c  if(x(npts).cq.alon.and.npts.gt.5)x(npts)  =  alonmx 

c  IF  (ZONE.GT.30)X(NPTS)=-XINTNUM/36000.  +  ALON 

lDX=lDX  +  36 
ENDDO 


wrilc(9,1604)ifcatp 

if(ifcatp.cq.0)thcn 

WRITE(9.1600)IFAC,IFICN,ISMC,IPHT. 

*  poricnt,plcngth,pwidth,NPTS 
cndif 

if(ifcatp.cq.l)thcn 

WRITE(9.1 60 1  )IFAC,IFICN,ISMC.IPHT. 

*  ldircct,lwidth,NPTS 
cndif 

if  (ifcatp.cq.2)thcn 

WRITE(9, 1 602)IFAC,IFICN.1SMC.IPHT. 

*  astruct,atrcc,aroof,NPTS 
cndif 

1600  format(8I5) 

1601  format(7I5) 

1602  formal(8I5) 

1604  fomiat(I2) 


DO  MM=1,NPTS 
WRITE(9,9876)X(MM).Y(MM) 
9876  FORMAT{2F15.10) 

ENDDO 


GO  TO  300 

Q  ilHI,^,it'*>K***************  **************************  ********* 

800  CONTINUE 

c  WRITE(9,3400)aHEAD-l),IREJECT,IPHTMAX 
3400  FORMATC  END  OF  MANUSCRIPT’;, 

*  -TOTAL  FACS  =’,I9,/, 

*  '  REJECTED  FACS  =',19;. 

*  '  MAXIMUM  FEATURE  HEIGHT  =’,19,/) 

C—  GET  NEXT  MANUSCRIPT 

c  idifr=600-i36+l 
c  print  ^.'Number  of  36  read  =  ',i36 
c  print  *, 'Reading  to  end  of  record  =  '.idiff 
c  do  i=l, idiff 
c  rcad(2,500)bits36 
c  enddo 

c  print  *,'Rcading  checksum  word' 
c  rcad{2,500)bits36 


GOTO  100 

C - 

C  ALL  MANUSCRIPTS  PROCESSED 
C - 


900  CONTINUE 
DO  1=1,1000 

IF  (lFEAIDa,l).GT.O  .OR.  IFEAID(L2).GT.0  .OR. 

*  IFEAID(I,3).GT.0) 

*  WRITE(6,2700)I,aFEAID(I,J),J=l,3) 
enddo 

WRITE(6, 1 700)0,(ISMCAR(I,J).J=1 ,3),I=1 , 1 4),1MISS 

CLOSE  (UN1T=2) 

CLOSE  {UNIT=3) 
close  (unit=9) 

STOP 

1 1 00  FORMATf//,’  MANUSCRIPT  NUMBER:', IS/ 

*  '  LEVEL  NUMBER  :',I5;. 

*  •  WAG(WAC)  NUMBER  :',I5;, 

*  •  WAG(WAC)  CELL 

*  •  WAG  CELL  :',I5;. 

*  •  SOUTHWEST  LON/LAT;',2FI  2.2;, 

*  '  NORTHEAST  LON/LAT:’,2F  1 2.2/) 

1200  FORMAT(/,’  DATA  SET  ID  RECORD: ',/, 

*  '  DATA  CLASSIFICATION  :’,Ai;, 

*  '  SECURITY  HANDLING  :',27A1,/. 

*  '  PRODUCT  TYPE  :'.5Ai;, 

*  '  MANUSCRIPT  REF.  NUMBER: ',15Ai;, 

*  ’  DATA  EDITION  NUMBER  : ',2Ai;. 

*  '  MATCH/MERGE  VERSION  :',Ai;, 

*  '  MAINTENANCE  DATE(YYMM): ',4 A  i;, 

*  '  MATCH/MERGE  DATE(YYMM):  ',4Ai;, 

*  '  PRODUCT  SPEC.  STOCK  NO: ',9A1,/. 

*  '  AMMENDMENT/CH ANGE  NO:  *,2 A 1 J, 

*  '  DATE(YYMM)  :'.4Ai;, 

*  *  HORIZONTAL  DATUM  CODE: ',5  A I ;, 

*  '  COMPILATION  DATE{YYMM):’.4A1,/) 

1300  FORMATC  ACCURACY  RECORD:’;, 

*4X,4A1,’  =  ABSOLUTE  HORIZONTAL  ACCURACY  (M)’;, 
MX,4A1,’  =  POINT-TO-POINT  HORIZONTAL  ACCURACY  (M)',/, 
*4X,4A1,'  =  VERTICAL  HEIGHTING  ACCURACY  (M)',/, 
*4X,2A1,'  =  MULTIPLE  ACCURACY  OUTLINE  FLAG';, 

*10X,’  00  =  NO  ACCURACY  SUBREGIONS’;, 

*10X, '02-09  =  NUMBER  OF  ACCURACY  SUBREGIONS';//) 


1400  FORMATflx,’ - Point - ';, 

*  2X,’FAC=',15,’  Fca  typc=',1 1 ,’  HT='.I4, 

*  2X,'FIC=’,I4,'  SMCC=’,I2,'  Num  coor  fca=’,I5, 

*  2X,’Ori=',I2,’  Len  =',I3,'  Width  =  ',13) 

1401  FORMAT(lx,' - Linear- - ’;, 
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*  2X,TAC=M5,'  Fca  typc=M  1 HT=M4, 

*  2X;FIC=M4,’  SMCC=M2;  Num  coor  fca=',15, 

*  2X;DiKn,'Widlh=M3) 

1402  FORMAT(lx; - Area - 7, 

*  2X;FAC=\I5;  Fca  typc-'.Il HT=M4, 

*  2X,'FIC=M4,'  SMCC=',I2.'  Num  coor  fca=M5, 

*  2X,'Strut=',I2;  Trcc=M3;  Roof^M3) 


1 700  FORMATC  SOIL  MATERIAL  TYPE  POINT  LINE  AREA7. 


1X,I5,' 

METAL 

’,315,/, 

1X,I5,’ 

PART  METAL  ',315,/, 

1X,I5.’ 

STONE 

',315/, 

IX, 15,’ 

COMPOSITION  ',3157 

1X,15,' 

EARTHWORKS  '.315,/, 

1X.15.’ 

WATER 

',315/, 

IX, 15.’ 

DESERT 

’.315/, 

IX, 15,' 

ROCK 

’,315,/, 

1X,15,’ 

CONCRETE  ’,315/, 

1X,15,' 

SOIL 

’,315/, 

IX, 15,' 

MARSH 

’.315/, 

1X,I5,’ 

TREES 

’,315/, 

1X,I5,' 

SNOW 

',315/, 

1X,15,’ 

ASPHALT 

’,315/, 

■»  •  NO  CODE  ',15) 

1800  FORMATC  NUMBER  OF  CELLS  WITH  VEG  HTS-M8) 
2700  FORMATC  FEATURE  NUMBERiPOINT  LINE  AREA'.4I5) 

END 


C-  12  =  number  ofbits  to  be  used 
C-  XINTNUM  =  NUMBER  RETURNED 


COMMON /A/ BUFF36{12,8),BUFFBlTS(288),buttbits(288) 
COMMON  /D/  ADATA(21 600*30) 

BYTE  BUFF36,BUFFBlTS,ADATA,buttbits 
INTEGER*4  Il,I2,I3,J,ICOUNT 
INTEGER*4  XINTNUM 

ICOUNT=-l 

XINTNUM=0 

IWORK=0 

13=12+11-1 

DO  J=13, 11,-1 
ICOUNT=ICOUNT+l 
ITEMP=ADATA(J) 

IF(1C0UNT.LT.32)CALL  MVBITSOTEMP.O,  1  .XINTNUM,ICOUNT) 
END  DO 

RETURN 

END 


SUBROUTINE  C32T36(BUFF) 

,i,t>****^********'t‘*  ********  *******<>‘*<1'*  ********** 

C- 

C-  This  subroutine  converts  nine  32  bit  words  into 
C-  eight  36  bit  words. 

C- 

COMMON  /A/  BUFF36(12,8).BUFFB1TS(288) 

BYTE  BUFF36,BUFFB1TS 
INTEGER*4  BUFF(9),AWORK,AOUT 


aout  =  0 
c  m=31 

DO  101=  1,288 
J  =  (M)/32  +  l 
K  =  I-I 

8  1F(K.LT.32)G0T0  9 
K  =  K  -  32 

GO  TO  8 

9  CONTINUE 
AWORK  =  BUFF(J) 

c  asavc  =  buffQ) 

1CK  =  3I  -K 

CALL  MVBITS(AWORK,KK.1,AOUT.O) 
BUFFBlTS(l)  =  AOUT 
c  call  mvbits(asavc,m,  1  ,aout,0) 

c  buttbits(i)  =  aout 

c  m=m-l 

c  if(m.cq.-l)m=31 

aout  =  0 

10  CONTINUE 
c  writc(9,100) 

c  writc(9,101)buff 
c  writc(9,110) 
c  writc(9,l  1  l)buffbits 
c  writc(9,120) 
c  writc(9,121)buttbits 


cl 00  format(lx,' - buff- - *) 

cl  10  format(Ix,’ - buffbits - *) 


cl 20  format(lx,’ - BUTTbits - ■) 

cl  01  format(9il6) 
cl  11  format(36i2) 
cl21  format(32i2) 

RETURN 

END 


SUBROUTINE  BlN2DEC_Dai, 12, XINTNUM) 


C- 

C-  Converts  to  decimal  (DELTA  REFERENCE) 
C-  11=  starting  bit  location 
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Appendix  G  -  PV-WAVE®  Procedures 


arics.com 

common  hdr,  dcg_to_rad,  rad_to_dcg,  alpha,  ccc_sq.  $ 
nc_origin,  rad„IIc_origin,xyz_origin_uvw,  $ 
gsmtxOl,  gsmtx02,  gsmtx21,  gsmtx22,  S 
rad_llc,uvw_conv,uvw_ofTsct,xyz_conv 


pro  makc_hdrbin,hfilc 

ilun  =  5 

gcLiun.ilun 

opcnr.ilun,hfilc 

f  =  fstat(ilun) 

h  =  bytarr(f.sizc) 

rcadu,ilun,h 

pointjun,ilun.0 

id9  =  whcrc(h  cq  9) 

id32  =  whcrc(h  cq  32) 

isz  =  strtrim(string(  id32(0)  -  id9(2)  - 1  ),2) 

nx^O 

ny  =  0 

fmtl  “'(19x,i'+  isz  +  ',lx,i'  +  isz  +  ',2x//)' 

rcadf,ilun,ny.nx,format=fmtl 

fmt2  =  ’(22x.i4,lx,f8.5,4x,i4.1x,R.5)' 

rcadf,iIun,swJat_d,sw_lat_m,sw_lon_d,sw_[on_m,format=fmt2 

rcadf,ilun,nc_lat_d,ncjat_m,nc_lon_d,ncjon_m.format=fmt2 

fmt3  ='(30x,f8.6,lx,f8.6) 

rcadf,ilun,row_scl,co!„scl,format=fmt3 

frcc_lun,iliin 

olun  =  5 
gct_lun.olun 
opcnw.olun,hfl  Ic+'.ncw' 
writcu.olun,nx,ny 

writcu.olun,sw_lat_d,sw_Iat_m,sw_ion_d,sw_lon_m 
write  u,olun,nc_lat_d,nc_Iat_m,nc_lon_d,nc_lon_m 
writcu,olun,row_scl,coLscI 
frcc_lun,olun 

return 

end 


pro  sct_valucs 
@arics.com 

radJlc_origin  =  {,1at:O.OD,  lon:O.OD,  clv:O.OD  ) 

dcg^to_rad  «  1 .745329251 99D-2 
rad_to_dcg  =  57.2957795 1 32D 
alpha  =  63781 37.0D 

ecc„sq  =  (2.0D -  (l.OD / 29 8. 25 7223 5 63 D))  *  (l.OD / 298. 257223563 D) 

rad_lIc_origin.Iat  =  llc_origin.lat  *  dcg_to_rdd 
radJlc_origin.lon  =  llc„origin.lon  *  dcg_to„nid 

sin_llc_originJat  =  sin(rad_llc_origin,lat) 
cos_llc„origin_lat  =  cos(rad_llc_origin.lat) 
sinJlc„origin_lon  =  sin(radJlc_origin.lon) 
cos_llc_originJon  =  cos(rad_llc_origin.!on) 

gsmtxOl  =  -sinJlc_origin_lat 
gsintx02  =  cosjlc_originjat 
gsmtx21  =  cosjlc_originjat 
gsmtx22  =  sin_lIc_origin_lat 

xyz_origin_uvw  =  llc2uvw(llc_origin) 

return 

end 


pro  tcs2binary 

rtyp  =  0 

nent  =  0 

xpt  =  0.0 

ypt  =  0.0 

cr  =  J5tring(''15b) 

xmin  =  999999.0 

ymin*  9999999.0 

xmax  =  -999999.0 


ymax  =  -9999999.0 

types  =  lonaiT(3) 

rfac  =  0 

rfic  =  0 

rsme  =  0 

rhgt  =  0 

rori  =  0 

rlcn  =  0 

rwid  =  0 

rpts  =  0 

rdir  =  0 

rstruc  =  0 

rtrcc  =  0 

rroof  =  0 

manu  =  1 

lastfac  =  -l 

fmtO  =  ’(8i5)'&szO  =  8*5+1 
fmtl  =’(7i5)'&szl  -7*5+1 
fmt2  =  '(8i5)'&sz2  =  8*5+1 
fmtdata  =  '(n3.3,n3.3,n3.3)'  &  szdata  =  13*3+1 

fnamc  = ' ' 

print, 'Change  TCS  Text  to  Binary’ 

print,' - - - ’ 

print.'Entcr  in  DFAD  file  to  process  -  NO  extension' 
rcad,fnamc 

opcnr.l , fnamc  +  *.tcs' 
opcnw,2,fname  +  '.tesbin' 

fin  =  fstat(l) 
fts  =  float(ftn.sizc) 

pent  =  0.0 

writcu,2,xmin,ymin,xmax,ymax 

while  not(cofl[l))do  begin 

rcadf,  1  .rtyp,format-'(i2)' 

writcu,2,rtyp 

pent  =  pent  +  2+1 

ease  rtyp  of 

0:  begin 

rcadf,  1 ,  rfac,rfic,rsmc,rhgt,rori,rlcn,rwid,ipts,format-fmtO 

writcu,2,rfac,rfic,rsmc,rhgt,rori,rlcn,rwid,rpts 

pent  =  pent  +  szO 

end 

1:  begin 

rcadfl ,  rfac,rfic,rsmc,rhgt,rdir,rwid,ipts,fonnat=frntl 

writcu,2,rfac,rfic,rsmc,rhgt,rdir,rwid,rpts 

pent  -  pent  +  szl 

end 

2:  begin 

rcadf,  I ,  rfac,rfic,rsmc,rhgt,rstruc,rtrcc,rroof,rpts,fonnat=fmt2 

writcu,2,rfac,rfic,rsmc,rhgt,rstruc,rtrcc,rroof,rpts 

pent  =  pent  +  sz2 

end 

endease 

txpts  =  fltarr(rpts) 
typts  =  fltarT(rpts) 
tzpts  =  fltarr(rpts) 
zpt  =  0.0 

for  i=0,rpts-l  do  begin 
rcadf,  1  .xpt, ypt,zpt,format-fmtdata 
txpts(i)  =  xpt 
typts(i)  =  ypt 
tzpts(i)  =  zpt 
pent  =  pent  +  szdata 
end 

types(rtyp)  =  typcs(rtyp)+(rtyp+l) 

writcu,2,txpts,typts,tzpts 

txmin  -  min(txpts,max-txmax) 
tymin  =  min(typts,max=tymax) 
if  txmin  It  xmin  then  xmin  =  txmin 
if  txmax  gt  xmax  then  xmax  -  txmax 
if  tymax  It  ymin  then  ymin  =  tymin 
if  tymax  gt  ymax  then  ymax  =  tymax 

if  rfac  It  lastfac  then  manu=manu+l 
lastfac  =  rfac 
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print, manu,rtyp.rfic,rpts,rhgt,  S 

((pcnt/fts)*100.0),cr,fonnat='(S,5i5,2x,f6.2,"%",a)' 

end 

print,' ' 

print, 'Number  of  point,  linear,  area  features’ 
print.typcs(0),typcs(l)/2,typcs(2)/3 

close,! 

point_lun,2,0 

writcu,2,xmin,ymin.xmax,ymax 

close, 2 

end 


pro  investigate 

rtyp  =  0 

nent  =  0 

xpt  =  0.0 

ypt  =  0.0 

cr=string("15b) 

xmin  =  999999.0 

ymin  =  9999999.0 

xmax  =  -999999.0 

ymax  =  -9999999.0 

types  =  lonaiT(3) 

rfac  =  0 

rfic  =  0 

rsme  =  0 

rhgt  =  0 

rori  =  0 

rlcn  =  0 

rwid  =  0 

rpts  =  0 

rdir  =  0 

rstruc  =  0 

rtrcc  =  0 

rroof  =  0 

manu  =  I 

lastfac  =  -1 

fnamc  = ' ' 


print,'Entcr  in  DFAD  file  to  investigate' 
rcad,fnamc 

opcnr.l , fnamc  +  '.dfadbin' 
opcnw,2,fnamc  +  '.check' 

fin=  fstal(l) 
fls  =  float(fin.sizc) 

pent  =  OL 
ekent  =  OL 
sbits  =  750L 

rcadu,  I  .xmin,ymin,xmax,yinax 

while  not(cof(l))do  begin 

rcadu,  Lrtyp 
pent  =  pent  +  2 
sbits  =  sbits  +  2 
ekent  =  ekent  +  2 

ease  rtyp  of 

0:  begin 

rcadu.1 ,  rfac,rfic,rsmc,rhgt,rori,rIcn,rwid,rpts 

pent  =  pent  +  (8‘*‘2) 

end 

1:  begin 

rcadu,  I ,  rfac,rfic,rsmc,rhgt,rdir,rwid,rpts 

pent  =  pent  +  (7‘*‘2) 

end 

2:  begin 

rcadu,  1 ,  rfac, rfic,rsmc,rhgt, rstruc, rtrcc,rroof, rpts 

pent  =  pent  +  (8'^2) 

end 

cndcasc 

txpts  =  fitantrpts) 
typts  =  fltanfrpts) 
rcadu,  l,txpts,typts 


pent  =  pent  +  (rpts  4  2) 
sbits  =  sbits  +  rpts 
ekent  =  ekent  +  rpts 

if  (ekent  gt  600)thcn  begin 
ekent  =  0 
sbits  =  sbits  +  2 
endif 

if  (rtyp  cq  0)thcn  goto,skip 

diffx  =  txpks(0:*)-txpts(l:'*) 

diffy  =  typts(0:*)-typts(l  :■") 

idx  =  whcrc(difTx  gt  1 .0  or  diffx  It  -1 .0.xcnt) 

idy  =  whcrc(diffy  gt  1 ,0  or  diffy  It  -1 .0.ycnt) 

if  (xent  gt  0  or  yent  gt  0)thcn  begin 
printf,2,'FAC,x,y,pts  ',rfac,xcnt,ycnt,rpts,sbits 
printf,2,'idx,idy  ',idx,idy 
if  (xent  gt  0)thcn  begin 

printf,2,'Xdiff  =  ',difrx(idx) 
printf,2,txpts(idx-l  ),typts(idx-l ) 
printf,2,’— ' 

printf,2,txpts(idx  ),typts(idx) 
printf,2,'— ' 

printf,2,txpts(idx+l),typts(idx+l) 

endif 

if  (yent  gt  0)thcn  printf,2,'Ydiff  =  ',difTy(idy) 

printf,2,' - ' 

endif 

skip: 

print,rfac,rtyp,rfic,rpts,  S 

((float(pcnt)/fts)’^  1 00,0),cr,format='(S,4i5,2x,f6.2,"%",a)' 


cndwhilc 

print,pcnt 
close,  1 
closc,2 

stop 

end 


pro  arics2binary 

rtyp  =  0 

nent  =  0 

xpt  =  0.0 

ypt  =  0.0 

cr=  string("15b) 

xmin  =  999999.0 

ymin  =  9999999.0 

xmax  =  -999999.0 

ymax  =  -9999999.0 

types  =  lonarr(3) 

rfac  -  0 

rfic  =  0 

rsmc  =  0 

rhgt  =  0 

rori  =  0 

rlcn  =  0 

rwid  =  0 

rpts  =  0 

rdir  =  0 

rstruc  =  0 

rtrcc  =  0 

rroof  =  0 

manu  =  1 

lastfac  =  -1 

fmtO  =  '(8i5)’  &  szO  =  8*5+1 
fmtl  =’(7i5)'&S2l  =7*5+1 
fint2  =  ’(8i5)'&sz2  =  8*5+1 
fmtdata  =  ’(fl5.10,fl5.10)’  &  szdata=  15*2+1 

fnamc  =  ’ ' 

print,'Changc  DFAD  Text  to  Binary' 

print,' - ’ 

print,'Entcr  in  DFAD  file  to  process  -  NO  extension' 
read, fnamc 

opcnr,l, fnamc  +  '.dfad' 
opcnw,2,fnamc  +  '.dfadbin' 

fin  =  fstat(l) 
fts  =  float(fin.sizc) 
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pent  =0.0 

wrilcu,2,xmin,ymin,xmax,ymax 

while  not(cof(l))do  begin 

rcadf,l  ,rtyp,format=’(i2)' 

writcu,2,rtyp 

pent  =  pent  +  2+1 

ease  rtyp  of 

0:  begin 

readf.l ,  rfac,rfic,rsmc,rhgt,rori,rlcn.rwid,ipts,fonnat=fmtO 

writcu,2,rfac,rfic,rsmc,rhgt.rori,rlen,rwid,rpts 

pent  =  pent  +  sz0 

end 

I;  begin 

rcadf,! ,  rfae,rfic,rsmc,rhgt,rdir,rwid,rpts,fomiat=fhitl 
writcu  ,2,rfae  ,rfic,rsmc,rhgt,rdir,rwid,rpts 
pent  =  pent  +  szl 
end 

2;  begin 

rcadf,) ,  rfac,rfic,rsmc,rhgt,rstruc,rtrce,rroof,rpts,format=fmt2 

writcu,2,rfac,rfic,rsme,rhgt,rstrue,rtrcc,rroof,rpts 

pent  =  pent  +  sz2 

end 

endease 

txpts  =  fltarr(rpts) 
typts  =  fltarr(rpts) 

for  i=0,rpts-I  do  begin 
readf.l  ,xpt,ypt, forma  t=fmtdata 
txpts(i)  =  xpt 
typts(i)  =  ypt 
pent  =  pent  +  szdata 
end 

typcs(rfyp)  =  typcs(rtyp)+(rtyp+l) 

writeu,2,txpts,typts 

txmin  =  min(txpte,max=txmax) 
tymin  =  min(typts,max=tymax) 
if  txmin  It  xmin  then  xmin  =  txmin 
if  txmax  gt  xmax  then  xmax  =  txmax 
if  tymax  It  ymin  then  ymin  =  tymin 
if  tymax  gt  ymax  then  ymax  =  tymax 

if  rfae  It  lastfac  then  manu=manu+l 
lastfae  =  rfae 

print, rfac,rtyp,rfic,ipts,rhgt.  S 

((pcnt/fts)*100.0),cr, formats, 5i5,2x,f6.2,”%", a)' 

end 

print,' ' 

print.'Number  of  point,  linear,  area  features' 
print ,  types  (0),types(  1  )/2,typcs(2)/3 

close,! 

point_lun,2,0 

writcu,2,xmin,ymin,xmax,ymax 

closc,2 

end 


;  Sid;  color_palcttc.pro,v  1,1  1991/05/22  16:53:19  jeffry  Exp  S 

pro  cpal 

;+ 

i  NAME:  COLOR.PALETTE 

;  PURPOSE:  To  display  the  numerical  values  associated  with  a  color  table 

;  Find  interval  to  be  used  on  the  table 

int=l 

inx=l 

if (!d.n_colors  gt  128)  then  inx=.5 

;  Save  currently  active  window  number 
holdw=  Id.window 
holdp  =  Ip.color 


;  Set  up  window 

yboxcs=fix(!d.n_colors/(8'''int)) 

yva!uc=yboxcs*8*int 

if  (yvalue  nc  !d.n_colors)  then  yboxcs=yboxcs+l 

ysizc=yboxcs'''40*inx 

window,frcc=  1  ,xsizc=320,ysizc=ysi2c 

;  Calculate  when  to  switch  printing  the  label  in  the  opposite  color 

!p.color=!d.n_colors-l 

changc=yvaluc'''3/4 

;  Loop  through  colors 

x=0 

y=0 

for  i=0,!d.n_colors-l,int  do  begin 
tv,rcplicatc(i,40,(40*inx)),x,y 
if(i  gc  change)  then  !p.color=0 
xyouts,x+5,y+5,strtrim(string(i), 2), /device 
x=x+40 

if  (x  ge  40’*8)  then  begin 
x=0 

y=y+(40*inx) 

endif 

endfor 

;  Set  back  to  previously  active  window 
if  holdw  gc  0  then  wsct.holdw 
!p.color=hoIdp 

return 

end 


pro  showbinary 

rtyp  =  0 

nent  =  0 

xpt  =  0.0 

ypt  =  0.0 

cr  =  string{"15b) 

xmin  =  999999.0 

ymin  =  9999999.0 

xmax  =  -999999.0 

ymax  =  -9999999.0 

types  =  lonarr(3) 

rfae  =  0 

rfic  =  0 

rsme  =  0 

rhgt  =  0 

rori  =  0 

rlcn  =  0 

rwid  =  0 

rpts  =  0 

rdir  =  0 

rstruc  =  0 

rtrcc  =  0 

rroof  =  0 

manu  =  1 

lastfac  =  -! 

fiesav  =  intarT(3,1000) 

ficcol  =  [0,0,2,6,0,7,16,23,8,3,30,  8,15, 4,27, 3] 
fnamc  = ' ' 

print.'Enter  in  DEAD  file  to  view* 
rcad.fnamc 

opcnr.l , fnamc  +  '.dfadbin' 
cfac  =  0 

print.'Entcr  FAC  to  view  or  zero  for  all' 
rcad.cfac 

pent  =  1 

rcadu,l,xmin,ymin,xmax,ymax 

dcvicc,pscudo_color=8 

window,0,xsizc=900,ysizc=900 


tck_color 

plot,[xmin,xmin,xmax,xmax,xmin],[ymin,ymax,ymax,ymin,ymin],xstylc=l,ystylc=l,$ 

color=0,tickformat='(f6.3)',background=l 

while  not(coftl))do  begin 

rcadu,l,rtyp 

ease  rtyp  of 
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0:  begin 

rcadu.l ,  rfac,rfic,rsmc,rhgt,rori,rlcn,rwid,rptj; 
end 

1 ;  begin 

rcadu,l ,  rfae,rfic,rsmc,rhgt,rdir,rwid,ipts 
end 

2:  begin 

rcadu,l,  rfac,rfic,rsmc,rhgt,rstruc,rtrcc,rroof,rpts 
end 

cndcasc 

txpts  =  fltarr(rpts) 
typts  =  fltarr(rpts) 
rcadu.l,  txpts,  typts 

if  rfac  nc  cfac  and  cfac  nc  0  then  goto.skipper 
ficidx  =  rfic/1 00 

if  (ficidx  cq  9)thcn  ficidx  =  10  +  ((rfic/10)-(rfic/100*  10)) 


radjic.lat  =  llc.Iat  *  dcg_to_rad 
radjle.lon  =  llc.lon  *  dcg_to„rad 

:  Convert  geodetic  to  JointStars-gcoccntric 


lon_ofrsct  =  radjle.lon  -  rad_llc_origin.lon; 
cos_lat  =  cos(  radjic.lat ); 
sinjat  =  sin(  radjic.lat ); 

re  =  alpha  /  sqrt(  1 .0  -  (ccc_sq  *  sinjat  *  sin_lat)  ) 

uvw_conv.x  =  (re  +  llc.clv)  *  cosjat  *  cos(  !on_offsct ) 
uv\v„conv.y  =  (re  +  llc.clv)  *  cosjat  *  sin(  !on_ofTsct ) 
uvw_conv.z  =  ( (re  *  (1 .0  -  ccc_sq))  +  llc.clv )  *  sinjat 

rctum,uvw_conv 

end 


if  rfac  cq  1  then  goto, skipper  - 

wid  =  rwid/2  function  uvw2tcs,uvw 

Icn  =  rlcn/2 

@arics.coni 

if  rtyp  cq  0  and  rfic  nc  420  then  oplot,txpts,typte,color=ficcol(ficidx),psym=l  ,symsizc=0.5 


if  rtyp  cq  0  and  rfic  gt  400  and  rfic  It  500  then  begin 
print, rwid.rlcn 

oplot,[txpts-wid,txpts-wid,txpts+wid,txpts-(-wid,txpts-wid],  $ 
[typts-lcn,typts+lcn,typts+lcn,typts-icn,typts-lcn],  S 
color=ficcol(ficidx) 
stop 

endif 

if  rtyp  cq  1  then  oplot,txpts,typts,colop=ficcol(ficidx) 
if  rtyp  cq  2  then  polyfill, txpts,typts,color=ficcol(ficidx) 
ficsav(rtyp,rfic)  =  ficsav(rtyp,rfic)  +  1 

if  rfac  cq  cfac  and  cfac  nc  0  then  begin 
stop 
endif 

pent  =  pent  +  1 

if  rfac  It  lastfac  then  begin 
manu=^manu+l 
print, manu 
endif 

lastfac  =  rfac 

skipper; 

end 

print, pent 
close,  1 


;  convert  uvw  to  tes 


uvw_ofTsct.x  =  uvw.x  -  xyz_origin_uvw.x 
uvw_offset.z  =  uvw.z  -  xyz_origin_uvw.z 

xyz_conv.x  =  uvw.y 

xyz_conv.y  =  gsmtxOl  *  uvw_offsct.x +  gsmtx21  *  uvw„ofrsct.z 
xyz_conv.z  =  gsmtx02  *  uvw_offsct,x  +  gsmtx22  *  uvw_offsct.z 

rctum,xyz_conv 

end 


pro  find_dfad_clcv 
@arics.com 

llc_origin  =  {,lat:0.0D,  lon;0.0D,  cIviO.OD  } 
lie  =  {,lat:0.0D,  lon:0.0D,  clv:0.0D  } 
xyz  =  {,  x:0.0D,  y:0.0D,  z:0.0D  } 
radjlc  =  {,lat:0.0D,  lon:0.0D,  clv:0,OD} 
uvw_conv  =  {,x:0.0D,  yrO.OD,  z:0.0D} 
uvw_offsct  =  {,x;0.OD,  y:0.0D,  z:0.OD) 
xyz_conv  =  {,x:0.0D,  y:0.0D,  z:0.0D} 


opcnw,2,fnainc  +  '.fie' 

printf,2.format=’(lx,"FIC’',5x,"Point",9x,"Uncar",8x,"Arca")’ 

printf,2,format=’(9x," - ",9x," - ",8x," — ")' 

formfic  =  "(Ix,i3,5x,i4,10x,i4,10x,i4)" 
for  i=l  00,999  do  begin 

if  ficsav(0,i)  nc  0  or  ficsav(l,i)  nc  0  or  ficsav(2,i)  nc  0  then  begin 
printf.2,i,ficsav(0,i),ficsav(l  ,i),ficsav(2,i),format=formfic 
endif 
end 

sumO  ==  long(total(ficsav(0,'^))) 
suml  « long(total(ficsav(l,*))) 
sum2  =  long(total(ficsav(2,*))) 
gtot  ==  sum0+suml+sum2 

printf,2,format='(lx,'Totals",2x," - ",9x," - ",8x," — ")' 

printf,2,suni0,suml,sum2,format='(9x,i5,9x,i5,9x,i5)' 
printf,2,gtot,format='(/,"Grand  total  =  ",i6)' 
close  ,2 


end 


function  llc2uvw,llc 

(garics.com 

Ion_offsct  -  O.OD 
cosjat  =  O.OD 
sin_lat  =  O.OD 
re  =  O.OD 


;  Decimal  degrees  to  radians 


print, 'Enter  in  Aries  file  to  find  elevations  for  -  NO  cxtcnsion.s' 
rcad.anamc 

openr,  1  ,anamc+'.dfad' 
opcnw,4,anamc+'.tcs' 

olat  =  29.000000D  &  olon  =  46.1 66666D  &  oclv  =  O.OD 

print, 'Enter  in  Origin  of  dataset  (Lat/Lon/Elcv)' 
llc„origin.lat  =  olat 
llc_origin.lon  ~  olon 
llc_origin.eIv  =  oclv 

; — Initialize  variables  for  LatLon  ->  TCS  conversion - 

sct_valucs 

edir  = '  ’ 

print, 'Enter  in  associated  elevation  directory' 
read, edir 
fnamc  = ' ' 

print, 'Enter  in  datafile  name' 
re  ad, fnamc 

hnamc  =  edir+'/hdr’ 
makc„hdrbin,hnamc 

openr, 2,cdir+'/'+fnamc 
openr, 3  ,edir+'/hdr.ncw' 


print, 'Reading  header  info . ’ 

xd  =  0&yd  =  0 
nx  =  0  «&;  ny  =  0 

swjat_d=0.0  &  swjatjm=0.0  &  swjon_d=0.0  &  swJon_m  =  0.0 
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ncJaLd=0.0  &  ncjat„m=0.0  &  nc_lon_d=0.0  &  ncJon„m  =  0.0 
row_scl=0.0  &  col_scl=0.0 


rcadu,3,nx,ny 

rcadu,3,swJat_d,swJaf_m,sw_lon_d.swJon_m,  S 
nc_lat_d,nc_lat_m,nc_lon_d,nc_lon_m,  S 
row_scI,coLscl 

sw_lon_scc  =  ( sw_lon_d*3600,0  +  sw_lon_m*60.0) 
sw_Iat_scc  =  ( sw_lat_d‘''3600.0  +s\v_lat_m*60.0) 
nc_lon_scc  =  (  nc_lon_d*3600.0  +  ncjon_m*60,0) 
nc_lat_scc  =  ( ncjat_d*3600.0  +  ncjat_m*60.0) 

print,'Crcating  elevation  array  NX  by  NY  ',nx.ny 
2  =  intacT(nx.ny) 

print, 'Reading  elevation  array . ‘ 

rcadu,2,z 


; - Processing  Text  file - ' 

rtyp  =  0  &  nent  =  0  &  xpt  =  0,0  &  ypt  =  0.0 
cr  =  string("15b) 

xmin  =  999999.0  &  ymin  =  xmin 
xmax  =  -999999.0  &  ymax  =  xmax 
types  =  lonarr(3) 

rfac  =  0  &  rfic  =  0  &  rsme  =  0  &  rhgt  =  0  &  rori  =  0  &  rlcn  =  0 
rwid  =  0  &  rpts  =  0  &  rdir  =  0  &  rstruc  =  0  &  rtrcc  ==  0  &  rroof  =  0 
manu  =  1  &  lastfac  =  -1  pent  =  0.0 

fmtO  =  '(8i5)’  &  szO  =  8*5+1 

fmtl  = '(715)’ &szl  =7*5+1 

fmt2  =  '(815)'  &  sz2  =  8*5+1 

fmtdata  =  '(HS.lO.flS.lO)'  &  szdata  =  15*2+1 

fmtout  ='(fl3.3,n3.3,n3.3)' 

fin  =  fsiat(l) 

Rs  =  floaKfm.sizc) 

while  not(cof(l))do  begin 

rcadf  ,1  ,rtyp,format='{i2)’ 
printf,4,rtyp,fomiat='(i2)' 
pent  =  pent  +  2+1 

ca.se  rtyp  of 

0:  begin 

rcadf  ,1 ,  rfac,rfic, rsme, rhgt, rori,rlcn,rwid,fpts,fomiat=fTntO 
printf,4,rfac,rfic,rsmc,rhgt,rori,rlcn,rwid,rpts,forTnat=fmtO 
pent  =  pent  +  szO 
end 

1 :  begin 

rcadf  ,1 ,  rfac,rfic,rsmc,rhgt,rdir,rwid,rpts,foiTnat=fmtl 
printf,4,rfac,rfic,rsmc,rhgt,rdir,rwid,rpts,format=fmtl 
pent  =  pent +  szl 
end 

2:  begin 

rcadf  ,1 ,  rfac,rfic,rsmc,rhgt,r<5truc,rtrcc,rToof,rpts,format=fmt2 
printf,4,rfac,rfic,rsmc,rhgt,rstruc,rtrcc,rroof,rpts,format=fmt2 
pent  =  pent  +  S22 
end 

cndcasc 

txpts  =  fltarr(rpts) 
typts  =  ntarr(rpts) 
tzpts  =  fltarr(rpt.s) 

for  i=0,rpt,s-1  do  begin 
rcadf,  1  ,xpt,ypt,format=fmtdata 
txpts(i)  =  xpt 
typts(i)  =  ypt 

xidx  =  fix(((xpt*3600.0)  -  swJon_scc)  /  col^scl) 
yidx  =  fix(((ypt*3600.0)  -  sw_lat_scc)  /  row_scl) 

if  (xpt*3600.0)  It  swJon_scc  or  $ 

(ypt*3600.0)  It  sw_lat_scc  or  $ 

(xpt*3600.0)  gt  nc_lon_scc  or  S 
(ypt  *3  600.0)  gt  ncjat_scc  then  begin 
tzptsfi)  =  z(0,0) 
endif  else  begin 

tzpts(i)  =z(xidx,yidx) 
cndclsc 

pent  =  pent  +  szdata 
end 


llc.lat  =  typts(i) 
llc.clv  =  tzpts(i) 
xyz  =  ]lc2tcs(llc) 
txpts(i)  =  xyz.x 
typts(i)  =  xyz.y 
tzpts(i)  =  xyz.z 
end 

for  i=0,ipts-l  do  begin 
xpt  =  txpts(i) 
ypt  =  typtsO) 
zpt  =  tzpts(i) 

printf,4,xpt,ypt,zpt,fonnat=fmtout 

end 

piint,rtyp,rfic,rpts.rhgt,((pcnt/fts)*100.0),cr,  $ 
format=’($,4i5,2x,f6.2,"%",a)' 


cndwhilc 

print,' ' 
close,  1 
closc,2 
closc,3 
closc,4 

end 


function  Ilc2tcs,lle 
@arics.com 

;  Convert  geodetic  to  JointStars-gcoccntric  (uvw) 
uvw  =  llc2uvw(llc) 

;  convert  JointStars-gcoccntric(uvw)  to  tes 

xyz  =  uvw2tcs(uvw) 

return, xyz 

end 


pro  showtes 

rtyp  =  0 

nent  =  0 

xpt  =  0.0 

ypt  =  0.0 

cr  =  string("15b) 

xmin  =  999999.0 

ymin  =  9999999.0 

xmax  =  -999999.0 

ymax  =  -9999999.0 

types  =  lonarr(3) 

rfac  =  0 

rfic  =  0 

rsme  =  0 

rhgt  =  0 

rori  =  0 

rlcn  =  0 

rwid  =  0 

rpts  =  0 

rdir=0 

rstruc  =  0 

rtrcc  =  0 

rroof  =  0 

manu  =  1 

la.stfac  =  -l 

ficsav=  intarr(3,1000) 

ficcol  =  [0.0,2.6.6,7,16,23,8,3,30,  8,15, 4,27,  3] 


fnamc  = ' ' 

print, 'Enter  in  TCS  file  to  view' 
read, fnamc 

opcnr.l , fnamc  +  '.tesbin' 
cfac  =  0 


for  i=0,rpts-l  do  begin 
lic.lon  =  txpts(i) 


pent  =  1 
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unzoom: 

rcadii,l,xmin,ymin,xmax,ymax 

zxmin  =  xmin 

zxmax  =  xmax 

zymin  =  ymin 

zymax  =  ymax 

Zoomit; 

dcvicc.pscudo_color=8 

window,0,xsizc=900,ysizc=900 


tck_color 

plot, [zxmin, zxmin, zxmax, zxmax,zxmin],  S 
[zymin,zymax,zymax,zymin, zymin],  S 
xstylc=l ,  ystylc=l ,  color=0,  tickformat^'(f6.3)',  background^! 

print,  zxmin, zymin, zxmax, zymax 

print.' - ' 

while  not(colTl  ))do  begin 

rcadu,l,rtyp 

ease  rtyp  of 

0;  begin 

ready,  l,rfac,rftc,rsmc,rhgt,rori,rlcn,rwid,rpts 
end 

1 :  begin 

rcadu,1 ,  rfac,rfic,rsmc,rhgt,rdir,rwid,rpts 
end 

2:  begin 

rcadu,l,  rfac,rfic,rsmc,rhgt,rstruc,rtrcc,rroof,rpts 
end 

cndcasc 

txpts  =  fltarr(rpts) 
typts  =  fltanfrpts) 
tzpts  =  fltarr(rpts) 
ready,  1 ,  txpts, typts,tzpts 


txmin  =  min(txpts,max=txmax) 
tymin  =  min(typts,max=tymax) 
if  txmin  It  zxmin  or  S 
tymin  It  zymin  or  S 
txmax  gt  kxmax  or  S 
tymax  gt  zymax  then  goto, skipper 


wid  =  rwid/2 
Icn  =  rlen/2 

if  rfac  nc  cfac  and  efac  nc  0  then  goto,skipper 
ficidx  =  rfic/lOO 

if  (ficidx  cq  9)then  ficidx  =  10  +  ({rfic/10)-(rfic/100*10)) 
if  rfac  cq  1  then  goto.skippcr 

oplot,txpts,typt.s,color=ficcol(ficidx),psym=l,symsi2c=0.5 
if  rtyp  cq  Olhen  begin 

oplot,[txpte-wid,txpts-wid,txpts+wid,txpts+wid,txpts-wid].  S 
[typts-Icn,typt.s+lcn,typts+Icn.typts-len,typts-lcn],  S 
color=ficcol(ficidx) 

endif 

if  rtyp  cq  1  then  oplot,txpts,typts,co!or=ficcol(ficidx) 
if  rtyp  cq  2  thcnpolyfill,txpts,typts,color=ficcol(ficidx) 
ficsav(iiyp,rfic)  =  ficsav(rtyp,rfic)  +  1 
if  rfac  cq  cfac  and  cfac  nc  0  then  begin 
stop 
endif 

pent  =  pent  +  1 

skipper: 

end 

print,pcnt 
ans  =  ” 

print.’Do  you  want  to  zoom  ?  y/n/all' 
rcad,ans 

if  ans  cq  'n'  then  goto,thccnd 
if  ans  cq  'all'  then  begin 
point_lun,l,0 


goto,unzoom 

endif 

print,'Selcct  first  point' 
cursor,xptl,yptiydata 
print,xptl,yptl 
wait.l 

print.'Sclcct  second  point' 

cursor,xpt2,ypt2ydata 

print,xpt2,ypt2 

zxmin  =  min([xptl  ,xpt2]) 

zymin  =  min((yptl  ,ypt2]) 

zxmax  =  max([xptl,xpt2]) 

zymax  =max([yptl,ypt2]) 

pointjun,l,0 

ready,  1  ,xmin,ymin,xmax, ymax 
goto,zoomit 

theend: 
close,  1 

opcnw,2,fnamc  +  '.fie' 

printf,2,fonnat='(lx,"FIC",5x,’'Point",9x, "Linear", 8x, "Area")' 

printf,2,format='(9x," - ",9x," - ",8x," — ")' 

formfic  =  "(Ix,i3,5x,i4,10x,i4,10x,i4)" 
for  1=100,999  do  begin 

if  ficsav(0,i) nc  0  or  ficsav(l,i)  nc  0  or  ficsav(2,i)  nc  0  then  begin 
printf,2,i,ficsav(0,i),ficsav(l,i),ficsav(2,i),format=formfic 
endif 
end 

sumO  =  long(total(ficsav(0,*))) 
suml  =  long{total(ficsav(l,*))) 
sum2  =  long(total(ficsav(2,*))) 
gtot  =  sum0+suml+sym2 

priritf2,format='(lx,’Totals",2x," - ",9x," - ",8x," — ")' 

printf,2,sum0,suml,sum2,format='(9x,i5.9x,i5,9x,i5)' 
printf.2,gtot,format='{/,"Gnind  total  =  ",i6)' 
close, 2 


end 


pro  dfad2tcs 
@arics.com 

llc„origin  =  {,!at;0.0D,  lon:0,0D,  clv:0.0D  } 

He"  =  {,lat:0.0D,  Ion;0.0D,  clv:0.0D  } 
xyz  =  {,  x:0.0D,  y:O.OD.  z:0.0D  } 
radjlc  =  {,lat;0,0D,  lon:0.0D,  clv;0.0D} 
uvw_conv  =  {,x:0.0D,  y:0.0D,  z:0.0D) 
uvw_offsct  =  {,x:0.0D,  y;0.0D,  2:0.00} 
xyz_conv  =  {,x:0.0D,  y:0.0D,  z:0.0D} 

anamc  = ' ' 

print,'Entcr  in  Aries  file  to  find  elevations  for  -  NO  extensions' 
read, anamc 

openr,  1  ,anamc+'.dfad' 
opcnw,4,anamc+'.tcs' 

olat  =  29.0000000  &  olon  =  46.1 66666D  &  ociv  =  0.00 

print, 'Enter  in  Origin  of  dataset  (Lat/Lon/Elcv)' 

llc_origin.!at  =  olat 
llc_origin.lon  =  olon 
llc_origin.elv  =  ocIv 

; — Initialize  variables  for  LatLon  ->  TCS  conversion - 

sct_valucs 

edir  = ' ' 

print,'Entcr  in  a.ssociatcd  elevation  directory' 

rcad,cdir 

fiiamc  = ' ' 

print,'Entcr  in  datafile  name' 
rcad,fnamc 

hnamc  =  edir+'/hdr' 
makc_hdrbin,hnamc 

opcnr,2,cdir+'/'+fnamc 
opcnr,3  ,cdir+'/hdr.ncw' 
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print.'Rcading  header  info . ' 

xd  =  0  &  yd  =  0 
nx  =  0  &  ny  =  0 

sw_lat_d=0.0  &  swjat_ni=0.0  &  sw_lon_d=0.0  &  swJon_m  =  0.0 
nc_lat_d=0.0  &  nc_lat_m=0.0  &  nc_lon_d~0.0  &  nc_Ion_rn  =  0.0 
row_scl=0.0  &  col_scl=0.0 


rcadu,3,nx,ny 

re adu,3  ,sw_l at_d,s w_! at_m ,s w_l on_d ,sw_! on_m,  S 
nc_lat_d,ncJat„ni,nc_lon_d,nc_lon_m,  S 
row_scl,col_scl 

sw_lon_scc  =  ( swJon_d*3  600.0  +  sw_lon_m*60.0) 
sw_lat_scc  =  ( swjat„d*3600.0  +  sw_lat_m*60.0) 
nc_lon_scc  =  ( nc_Ion„d*3600.0  +  nc_lon_tn*60.0) 
nc_lat_scc  =  (  ncjat_d*3600.0  +  ncjat_m*60.0) 

print.'Crcating  elevation  array  NX  by  NY  ',nx,ny 
z  =  intarr(nx,ny) 

print  .’Reading  elevation  array . ’ 

rcadu,2,z 

; - Proccs-sing  Text  file - ' 

rtyp  =  0  &  nent  =  0  &  xpt  =  0.0  &  ypt  =  0,0 
cr  =  string("15b) 

xmin  =  999999.0  &  yniin  =  xmin 
xmax  =  -999999.0  &  ymax  =  xmax 
types  =  lonarr(3) 

rfac  =  0  &  rfic  =  0  &  rsme  =  0  «Sr  rhgt  =  0  &  rori  -  0  &  ricn  =  0 
rwid  =  0  &  rpts  =  0  rdir  =  0  &  rstruc  =  0  &  rtrcc  =  0  &  rroof  =  0 
manu  =  1  &  lastfac  =  -1  &  pent  =  0.0 

fmt0  =  '(8i5)'&sz0  =  8*5+1 

fmtl  =’(7i5)'&szl  =  7*5+1 

fmt2  =  '(SiS)’  &  sz2  =  8*5+1 

fmtdata  =  '(fl5.10,n5.10)'  &  szdata  =  15*2+1 

fmtout  =’(n3.3,n3.3,n3.3)’ 

fin  =  fstat(]) 

fts  =  float(fin.size) 

while  not(cof(I  ))do  begin 

rcadf  ,1  ,rtyp,format='(i2)’ 
print  f,4,rtyp,format='(i2)' 
pent  =  pent +  2+1 

ease  rtyp  of 

0:  begin 

rcadf  ,1 ,  rfae.rfic, rsme, rhgt, rori, rlcn, rwid, ipts,format=fmt0 
printf,4,rfac,rfic,rsmc, rhgt, rori, rlcn,rvvid,rpts,fomiat=fnitO 
pent  =  pent  +  szO 
end 

1 ;  begin 

rcadf  .1 ,  rfac,rfic,rsmc,rhgt,rdir,r\vid,rpts,fonnat=fnitl 
printf,4,rfac,rfic,rsinc,rhgt,rdir,rwid,rpts,format=fmtl 
pent  =  pent  +  szl 
end 

2:  begin 

rcadf ,  1 ,  rfac,rfic,rsmc,rhgt, rstruc, rtrcc,rroof,rpts,format=fint2 
printf,4,rfac.rfic,rsmc,rhgt,rstruc,rtrcc,rroof,rpts,format=fmt2 
pent  =  pent  +  sz2 
end 

cndcasc 

txpts  =  ntarr(rpts) 
typts  =  fltaiT(rpts) 
tzpts  fltarr{rpts) 

for  i=0,rpts-l  do  begin 
rcadf,  1  ,xpt,ypt,fonnat=fmidata 
txpts(i)  =  xpt 
typts(i)  =  ypt 

xidx  =  fix(((xpt*3600.0)  -  sw_lon_scc)  /  coLscl) 
yidx  =  fix(((ypt*3600.0)  -  swjat^scc)  /  row_scl) 

if  (xpt*3600.0)  It  sw_lon_scc  or  $ 

(ypt*3600.0)  It  sw_lat_scc  or  $ 

(xpt*3600,0)  gt  nc_lon_scc  or  $ 

(ypt*3600.0)  gt  ncjat_scc  then  begin 
tzpts(i)  =  z(0,0) 
endif  else  begin 

tzpts(t)  =  z(xidx,yidx) 
cndcisc 

pent  =  pent  +  szdata 


end 


for  i=0,rpts-l  do  begin 
lie.  Ion  =  txpts(i) 
llc.lat  =  typts(i) 
llc.clv  =  tzpts(i) 
xyz  =  llc2tcs(lle) 
txpts(i)  =  xyz.x 
typts(i)  =  xyz.y 
tzptsO)  =  xyz.z 
end 

for  i=0,ipts-l  do  begin 
xpt  =  txpte(i) 
ypt  =  typts(i) 

2pt  =  tzptsfi) 

printf,4,xpt,ypt,zpt,format=fmtout 

end 

print,rtyp,rfic,rpts,rhgt,((pcnt/fts)*  1 00.0),cr,  S 
forTnat=’(S,4i5,2x,f6.2,'’%”,a)’ 


endwhilc 

print,'  ’ 
close,  1 
close, 2 
closc,3 
close, 4 

end 


pro  showall 

rtyp  =  0 

nent  ==  0 

xpt  =  0.0 

ypt  -  0.0 

cr  =  string("15b) 

xmin  =  999999.0 

ymin  =  9999999.0 

xmax  =  -999999.0 

ymax  =  -9999999.0 

types  =  lonarr(3) 

rfac  =  0 

rfic  =  0 

rsme  =  0 

rhgt  =  0 

rori  =  0 

rlcn  =  0 

rwid  =  0 

rpts  =  0 

rdir  =  0 

rstruc  =  0 

rtrcc  =  0 

rroof  =  0 

manu  =  1 

lastfac  =  -1 

fiesav  =  intarr(3,1000) 

ficcol  =  (0,0,2,6.6,7,16,23,8,3,30, 8,15, 4.27, 3] 


fnamc  =  ['aricsr,'arics2','arics3','arics4’,’arics5'] 


dcvicc,pscudo_color=8 

window,0,xsizc=900,ysizc=900 

xmin  =  43.000 
ymin  =  26.000 
xmax  =  49.000 
ymax  =  32.000 

tck_coIor 

plot,[xmin,xmin,xmax,xmax,xmin],[ymin,ymax,ymax,ymin,ymin],xstylc=l,ystylc=l,S 

color=0,tickformat='(f6.3)',background=l 


for  afilc  =  0,4  do  begin 


opcnr.l,fhamc(afile)  +  '.dfadbin' 
cfac  =  0 
pent  =  1 

rcadu.l  ,xmin,ymin,xmax,ymax 
print,fhamc(afilc),xmin,ymin,xmax,ymax 
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while  not(cof(l))do  begin 


rcadu,l,rtyp 
ease  rtyp  of 
0:  begin 

rcadu.l ,  rfac,rfic,rsmc,rhgt,rori,rlcn,rwid,rpts 
end 

1:  begin 

rcadu.l,  rfac.rfic.rsmc.rhgt.rdir.rwid.rpts 
end 

2:  begin 

rcadu.l ,  rfac.rfic.rsmc.rhgt.rstruc.rtrcc.rroof.rpts 
end 

cndcasc 

txpts  =  fltarr(rpts) 
typts  =  fltarr(rpts) 
rcadu.l, txpts,  typts 

if  rfac  nc  cfac  and  cfac  nc  0  then  goto,skippcr 
ficidx  =  rfic/100 

if  (ficidx  cq  9)thcn  ficidx  =  10  +  ((rfic/10)-(rfic/!  00*10)) 
if  rfac  cq  1  then  goto,skippcr 

if  rtyp  cq  0  then  oplot,txpts,typts,color=ficcol(ficidx),psym=l,synisizc=0,5 
if  rtyp  cq  1  then  oplot,txpts,typts,color=ficcol(ficidx) 
if  rtyp  cq  2  then  polyfill.txpts, typts, coIot=ficcol(ficidx) 
ficsav(rtyp,rfic)  =  ficsav(rtyp,rfic)  +  1 

if  rfac  cq  cfac  and  cfac  nc  0  then  begin 
stop 
endif 

pent  =  pent  +  1 

lastfac  =  rfac 

skipper: 

end 

close,  I 

olun=  0 
gcMun.olun 

openw,olun,fnamc(afile)+'.fic' 

printf.olun, form  at='(lx,'’FIC'’,5x, "Point", 9x,"Lincar",8x,"  Area")' 

printf,olun,fonnat=’(9x." - ",9x," - ".Bx," — ")' 

formfic  =  "(Ix,i3,5x,i4,10x,i4,10x,i4)" 
for  i=l  00,999  do  begin 

if  ficsav(0,i)  nc  0  or  ficsav(l,i)  nc  0  or  ficsav(2,i)  nc  0  then  begin 
printf,olun,i,ficsav(0,i),ficsav(l,i),ficsav(2,i),format=fomific 
endif 
end 

sumO  =  long(total(ficsav(0,*))) 
suml  =  )ong(total(ficsav(l,*))) 
sum2  =  long(total(ficsav(2,*))) 
glot  =  sumO+suml  +sum2 

printf,olun,fonnat='(lx,"TotaIs",2x," - ",9x,'' - '',8x," — ")' 

prinlf.olun,sum0,suml,sum2,format='(9x,i5,9x,i5,9x,i5)' 
printf.olun,gtot,fomiat='(/," Grand  total  =  ",i6)’ 
frccjun.olun 
ficsav(*,*)  =  0 

endfor 


end 


Appendix  H  -  Glossary 


2- D 

3- D 
ADS 
AFATDS 
AML 
API 

ARC/INFO® 


ArcView® 

ARIES 

AUTOGRAPHICS® 

C 

C^ISR 

CD-ROM 

DFAD 

DIS 

DT&E 

DIED 

ETE 

Excel 

FID 

FORTRAN 

GIS 

GUI 

ICD 

ID 

IEEE 

ERDAS  IMAGINE® 

JADS 

Janus 

Joint  STARS 

JTF 

LGSM 

LMTDS 

MIL-PRF 

MIL-STD 


2- dimensional 

3 - dimensional 

advanced  distributed  simulation 
Advanced  Field  Artillery  Tactical  Data  System 
ARC/INFO®  Macro  Language 
application  programmer  interface 

a  workstation  GIS  software  package  by  Environmental  Systems  Research 
Institute 

a  desktop  GIS  software  package  by  Environmental  Systems  Research 
Institute 

Advanced  Radar  Imaging  Emulation  System  developed  by  Lockheed 
Martin  Tactical  Defense  Systems,  Litchfield  Park,  Arizona 
a  GIS  software  package  by  Lockheed  Martin  Tactical  Defense  Systems, 
Akron,  Ohio 

a  coding  system  for  programming  scientific  problems  to  be  solved  by  a 
computer 

command,  control,  communications,  computers,  intelligence,  surveillance 
and  reconnaissance 

compact  disk  that  can  hold  a  large  quantity  of  computer  data 

digital  feature  analysis  data 

distributed  interactive  simulation 

developmental  test  and  evaluation 

digital  terrain  elevation  data 

End-To-End 

a  spreadsheet  application  by  Microsoft® 
feature  identification 

a  coding  system  for  programming  scientific  problems  to  be  solved  by  a 
computer 

geographic  information  system 
graphical  user  interface 
interface  control  document 
identification 

Institute  of  Electrical  and  Electronics  Engineers 

a  geographic  imaging  suite  by  ERDAS®,  Incorporated 

Joint  Advanced  Distributed  Simulation,  Albuquerque,  New  Mexico 

interactive,  computer-based  simulation  of  combat  operations 

Joint  Surveillance  Target  Attack  Radar  System 

joint  test  force  or  Joint  Test  Force,  Albuquerque,  New  Mexico 

light  ground  station  module 

Lockheed  Martin  Tactical  Defense  Systems 

military  performance  specification 

military  standard 
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MTI 

NIMA 

OT&E 

PDU 

PV-WAVE® 

R2V™ 

RPS 

RWS 

SAR 

SEDRIS 

STARS 

SWA 

T«&E 

TAFSM 

TCS 

TEXCOM 

TRAC 

TRADOC 

UTM 

VPF 

VSTARS 

WSMR 

WWW 


moving  target  indicator 
National  Imagery  and  Mapping  Agency 
operational  test  and  evaluation 
protocol  data  unit 

a  visual  data  analysis  software  package  by  Visual  Numerics,  Incorporated 
a  raster  to  vector  conversion  software  package  by  Able  Software 
Company 

radar  processor  simulation  developed  by  Northrop  Grumman,  Melbourne, 
Florida 

remote  workstation 
synthetic  aperture  radar 

Synthetic  Environment  Data  Representation  &  Interchange  Specification 

surveillance  target  attack  radar  system 

Southwest  Asia 

test  and  evaluation 

Tactical  Army  Fire  Support  Model 

Topocentric  Coordinate  System 

U.S.  Army  Test  and  Experimentation  Command 

U.S.  Army  Training  and  Doctrine  Command  (TRADOC)  Analysis  Center 

U.S.  Army  Training  and  Doctrine  Command 

Universal  Transverse  Mercator 

Vector  Product  Format 

Virtual  Surveillance  Target  Attack  Radar  System 
White  Sands  Missile  Range,  New  Mexico 
world  wide  web 
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deg  or  ° 

GB 

MB 


degree 

gigabyte 

megabyte 


Units  of  Measure 
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